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Abstract 

^ ■ The study of Finite Size Effects in Quantum Field Theory allows the extraction 

. of precious perturbative and non-perturbative information. The use of scaling func- 

r\ \ tions can connect the particle content (scattering theory formulation) of a QFT to its 

' ultraviolet Conformal Field Theory content. If the model is integrable, a method of 

investigation through a nonlinear integral equation equivalent to Bethe Ansatz and 
deducible from a light-cone lattice regularization is available. It allows to reconstruct 
the S-matrix and to understand the locality properties in terms of Bethe root configu- 
rations, thanks to the link to ultraviolet CFT guaranteed by the exact determination 
of scaling function. This method is illustrated in practice for Sine-Gordon / massive 
^ ■ Thirring models, clarifying their locality structure and the issues of equivalence be- 

tween the two models. By restriction of the Sine-Gordon model it is also possible to 
control the scaling functions of minimal models perturbed by ^(1,3). 
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Finite Size Effects (FSE) are widely recognized to play an important role in modern Sta- 
tistical Mechanics and Quantum Field Theory (QFT). From a statistical point of view, 
it is known that no phase transitions take place in a finite volume system. For example, 
the specific heat c(T), that is divergent at the critical point in infinite volume, looses it 
divergence if the system has finite size; one observes only a rounded peak in the plot of 
c(T) versus T (see fig.||). Moreover, there is only an interval around critical temperature 
Tc where the FSE are relevant. Out of this interval, they are negligible (because only near 
Tc the correlation length is comparable with the size of the system). The interesting fact 
is that specific heat (and other critical quantities) have a scaling behaviour (i.e. varying 
the size L) that is fixed by the (infinite size) critical exponents (see ||l|). This is a general 
fact: as argued in ||2|, the UV behaviour of the scaling functions (see later) is fixed by the 
conformal dimensions of the operators that belong to the universality class of the critical 
point (i.e. the CFT describing the critical point of the statistical system). 

Also in QFT interesting phenomena appear. If the space-time geometry is a cylinder of 
circumference L, Casimir effects change the energy of a two body interaction, because par- 
ticles interact in the two possible directions, as shown in fig. ||. New radiative corrections 
to the self-energy of a propagating particle may appear because of the closed geometry. 
For example, a massive particle propagating in a finite 1+1 dimensional space-time, with 
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Figure 1: Specific Heat near a phase transition point in infinite volume (dashed hue) and in finite 
volume (solid line) 

periodic boundary conditions in the space direction and an infinitely extended time direc- 
tion will have its propagator renormalized by the contributions of radiative corrections, 
as depicted in the second cylinder of fig.||. While on an infinite space only local virtual 
emissions are allowed, on a cylindrical geometry one can conceive virtual emissions that 
travel around the whole cylinder before coming back to the bare particle. Of course they 
are exponentially depressed if compared to the traditional local emissions, therefore, if 
the cylinder is large, their contributions to the self-energy of the propagator is negligible. 
However, if the cylinder circumference L is small enough, there is room left by Heisenberg 
principle for such virtual object to travel around the world and come back from the other 
side, thus giving new contributions that modify the self-energy. 

Liischer [|, || has estimated the relevance of such a phenomenon by computing the 
corrections to the mass of the particle due to FSE. For a massive self-interacting boson 
with (j)^ potential in 3+1 dimensions, they are represented by 

AM(L) = ^— [ dye-^'^^+y^Fiiy) + 0(e-v^™^) 

IGtt^ mL J 

m is the physical mass of the particle in infinite volume; F is the analytic continuation of 
the forward elastic scattering amplitude S with incoming momenta P11P2 and outcoming 
momenta P3,P4 

F{u) = S{pi,P2\p3,P4) 

where u = ^^p^^'^^p^^'P^'P^ and uj{p) = sjrr? ^ p^ . Saddle point integration peaks a contri- 
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bution from around y = 0, so that the result is 

AM oc — ^= + ... 

A very useful application of this estimate shows up in lattice QFT calculations. For 
example consider a 4 dimensional lattice, of lattice spacing a, on which you want to perform 
Montecarlo simulations of QCD. The lattice sizes are T x L x L x L, with, say, L = 20a. In 
lattice QFT, the discretization of space is to be intended as a regulator, i.e. it is equivalent 
to the introduction of a momentum cutoff A^ut ~ ^ which must be much larger than the 
typical hadron masses, if we want to get reliable data from the lattice. For example, we 
can choose a < 0.1 fm (1 fm = 10~^^ cm); then L = 20a = 2 fm, which is comparable with 
the proton radius (Compton wavelength). The proton in such a lattice is "squeezed", which 
means that the FSE are relevant. One way to overcome this difficulty is to build larger and 
larger computers, which, as history of the last 20 years indicates, goes to introduce other 
sorts of difficulties of financial, technical and also principal nature. The other possibility, 
at the cost of a pencil, is to be able to compute the FSE scaling of quantities related to the 
proton, from which we can extract physical data even from such small lattices. Of course 
nobody is nowadays able to estimate with accuracy the L dependence of quantities in a 
formidable problem like QCD, but progresses have been made in two dimensional QFT's, 
some of them sharing interesting phenomena with QCD or other interesting realistic 4 
dimensional QFT's. In particular, there are methods to exactly calculate the L dependence 
of physical quantities in the class of Integrable two dimensional QFT. This will be the main 
topic of the present review. 

Let us continue to analyze this lattice QFT example. Finite Size with periodic boundary 
conditions implies quantization of momenta in units of order Ap = ^ > 600 MeV. Pions 
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in such a lattice (m^r = 139 MeV) have virtual processes strongly affected by FSE. This 
illustrates a bit more numerically the observation about virtual processes made above. 
This discussion leads us to some important observations: 

• As momenta are quantized, the energy spectrum for finite L is discrete 

• FSE do not depend strongly on the lattice spacing a, they depend on the scale L ^ a 
and become irrelevant when L ^ Compton wavelength of the lightest particle (of 
mass m). Therefore the physical quantities may be expressed, up to overall factors 
encoding their dimensionality, as dimensionless functions of the dimensionless scale 
/ = mL. To such objects the usual denomination of scaling functions is given. 

• As L ^ a, FSE are relevant to renormalized QFT. / can play the role of dimensionless 
scale parameter of a sort of Renormalization Group fiow as we go to illustrate by the 
following heuristic arguments: 

— the limit / — > means or L i.e. the cylinder shrinks to a line, and there is 
no spatial dimension anymore, which is an uninteresting unphysical situation, 
or m — 0, which means that the lightest particle mass becomes irrelevant 
compared to the energies in play in the system, which is exactly the setup of the 
ultraviolet (UV) regime. Therefore, the limit / reproduces the UV regime 
of the model. 

— the limit / oo means or m — > oo, i.e. the situation where the masses of 
the particles are much greater than the scattering transferred momenta, which 
defines the infrared (IR) region, or L — ^ oo, i.e. the reconstruction of the QFT on 
the infinite volume. Note that these two limits are not in contradiction: actually 
the theory on the whole infinite volume is described in terms of the S-matrix, 
which is an intrinsically IR object, in the sense that it describes the scattering 
of asymptotic particles without looking at the details of the interaction at short 
range. 

• Both the UV and IR limits thus identified are fixed points of the QFT where the 
Callan - Symanzyk beta function annullates. They are therefore scale invariant 
points. Scale invariance in relativistic QFT often implies conformal invariance. 
Therefore the UV and IR points of the flow obtained by varying I from to oo 
are to be identified with Conformal Field Theories (CFT). 

— In the case of massive QFT, the IR CFT is trivial. The CFT description is useful 
in this case only at the UV point, where the theory can be seen as a relevant 
perturbation of a given CFT. At IR instead, the most effective description is in 
terms of the S-matrix. 

— In the case of massless (but not scale invariant) theories, both the UV and IR 
points are nontrivial CFT's, and the scaling functions interpolate between the 
physical quantities of these two CFT's. 
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2 Excursus on the Renormalization Group in 1+1 di- 
mensions 



In this section we quickly recall the main properties of a two-dimensional QFT under scale 
transformations. In general a QFT is not invariant under scale transformations, due to the 
presence of dimensionful parameters, appearing already at the classical level or generated 
by the renormalization procedure. To investigate this phenomenon it is useful to think in 
terms of the so called (Wilson) Space of Actions 0, |^. Let us suppose that the action 
S = f d^xC of a theory depends on a set of fundamental fields and their derivatives, 
and on a set of dimensionless coupling constants g = {gi,...,gn). The variation of the 
Lagrangian density C{gi, ■■■,gn) under scale transformations x ^ x + dtx can be seen as 
a transformation of the coupling constants g ^ g'. Scale invariance in this context means 
that under a scale transformation the point g remains unchanged. Scale transformations 
induce a sort of motion in the n-dimensional space of couplings. Such a space is called 
usually the Space of Actions. A trajectory in this space is a function g(t) (running coupling 
constants). The variation of the Lagrangian is C{t + dt) = C{t) + d^J^dt where = x^T^'-' 
is the Nother current associated to the scale transformation. The object d^J^ = = 6 is 
the trace of the stress energy tensor (we shall always consider in the following a symmetrized 
improved stress energy tensor 0) and scale invariance, where satisfies a continuity 
equation 9^ J'^ = 0, implies 9 = 0. For a general trajectory g(t) the stress energy tensor is 
not traceless. The Callan - Symanzyk beta function is defined as 

m ^ t (1) 

and can be interpreted as a sort of velocity field of the theory in the action space. Let 
us consider one coupling for simplicity. Imagine to be at a certain point in the space of 
actions where the action of the theory is characterized by a parameter g^. Consider this 
as a sort of initial condition, and study the behaviour of the theory as t increases. 

• if l3{go) > then g{t) is increasing in a region near go; 

• if P{go) < then g{t) is decreasing in a region near go. 

Points g* where P{g*) = are called fixed points. They divide the space of actions (one- 
dimensional in this example) in distinct regions of the g parameter. The running coupling 
constant keeps confined in one of these regions. Independently of the exact value of go, 
in a given region, the asymptotic value of g{t) for t — > ±oo is always given by the closest 
values of g where /3{g) = 0. There are two types of fixed points: 

• a fixed point g* is infrared (IR) if, starting from a value go of the coupling near it, 
one reaches it for t — > +oo. {P'{g*) < 0). 

• a fixed point g* is ultraviolet (UV) if, starting from a value go of the coupling near 
it, one reaches it for t — oo. {P'{g*) > 0). 
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In the case with more parameters all the beta functions have to annullate simultaneously 
in a fixed point. Given the initial conditions g = go there exists a unique trajectory, as 
the equation (|l|) is a first order linear differential equation. The trajectory g{t) is often 
referred in the literature as the Renormalization Group (RG) flow. 

The knowledge of the beta function is enough to reconstruct the behaviour of the theory 
in a surrounding of the point g^ i.e. to have a complete description of the tangent space 
to the space of actions in the neighborhood of g^. In such tangent space, the variation of 
the Lagrangian can be expressed as a linear combination of a base of fields identified as 
the derivatives of the Lagrangian with respect to the coupling constants gi 

ogi 

The trace of the stress energy tensor is a field living in such tangent space 

as it is a linear combination of the fields with coefficients given by the components of 
the beta function vector. It is then obvious that 6 = when and only when Pi{g) = 0, 
which means that the fixed points of a QFT are scale invariant. 

We go now to get an equation describing the variation of the point correlation 
functions of generic fields Ai{x) of the theory (X) = {Ai{xi), An{xn)) along the RG 
fiow. It is known in the literature as the Callan-Symanzyk equation. We present here a 
derivation inspired by A. Zamolodchikov The variation of correlation functions under 
scale transformation can be described in two equivalent ways 

1. apply the variation to all fields in the correlator 

6{X) = 6{A,{xi),...,AN{xM))=SjV^A,...A^e-^^^^ 
= 6JV^ [J2kAi...6Ak...ANe-^ -Ai...ANe-^6S] 

The variation of a generic field A^ under scale transformation is given by the global 
dilation of coordinates and by the intrinsic variation of the field due to its ingeneering 
(classical) dimension 

SAk = dtix'idl + Dk)Au 

The variation of the action is actually given by the insertion of the trace of the stress 
energy tensor into the correlator 

5S = j d'yQiy) 

Therefore 

S{X) = dt + Dk){X) - J d'y{Xe{y)) 
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2. alternatively one can vary the fields and the action with respect to the couplings gi 



but 

Be careful, in these formulae, not to confuse summations over i, that run in the 
tangent to the space of actions and summations over k that run over the points Xk 
in the correlator. Putting together the pieces, one gets 

i k '' 

Combining the two methods one arrives at the so called Callan - Symanzyk equation for 
the 1+1 dimensional RG of QFT 



(X) = 



where TkAk{xk) = [Ok + Xli Ak{xk) is the so called anomalous dimension of the 

field Ak. 

In particular the dimensions of the perturbing fields $i are given by 



|-(2L + 0) = E, (25^-S)$,(x) 



Applying this to the trace of the stress energy tensor 6 = S^^s r6 = 26. 

This, combined with Lorentz covariance, implies that also the other components of the 
stress energy tensor have the same anomalous dimension 2. This is also the ingeneering 
dimension of the stress energy tensor, which means that the stress energy tensor does not 
develop an anomalous dimension and keeps its classical dimension also at the quantum 
level. This is a fact common to all conserved currents. 

The most exciting result known in the study of two-dimensional Renormalization Group 
is undoubtedly the Zamolodchikov c-theorem 0. Here we do not reproduce the entire 
argument leading to this theorem (see the original paper § for a demonstration). We only 
offer the statement of the theorem in its full generality. 

c-theorem: In any unitary QFT2 a function c{g) exists, such that: 

1. it is decreasing along the RG flow: ^ < 0; 
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2. it is stationary in the fixed points g = g*, i.e. ^|^^ = if and only if ^ is a fixed 
point of the RG flow, i.e. Pi{g*) = 0. 

3. in such fixed points g*, the function c{g) fixes completely the two-point correlators 
of the stress energy tensor. Introducing the coordinates 

Z = Xi + 1X2 Z = Xi — 1X2 

and the notation T = T^z, T = Tgz, © = T^z = Tzz = T^^ one has = because the 
fixed point is scale invariant and 

(T(^, ^)T(0, 0)) = ^ (f (z, z-)f (0, 0)) = ^ (Tf ) = (2) 

where c = Yimfj^^ c{g) . The constant c is often called central charge and is a char- 
acteristic of the theory at the fixed point. 

As we have seen that theories are RG fiows from UV to IR fixed points when t goes from 
— oo to +00. This theorem means that, for unitary theories 

cuv > Cm 

which is the "simplified" form frequently mentioned in the literature. 



3 Conformal Field Theory 

In the previous section we have seen that the RG fixed points correspond to values of 
the parameters of a theory where scale invariance is guaranteed. Scale invariance implies 
under quite general assumptions a larger invariance, known as conformal invariance. In 
this section we summarize the main results in the treatment of conformal invariant theories 
(CFT) in two-dimensions, as developed in the fundamental work [in . 



3.1 The Virasoro algebra and the Hilbert space 

In two dimensions any analytic transformation of co-ordinates z = Xi+ix2 and z = Xi — ix2 
is a conformal transformation. The improved stress energy tensor obeys the continuity 
equation d^T^^ = and is symmetric T^'^ = T'^^. If we add the requirement of scale 
invariance, i.e. of null trace T^'" = 6 = 0, we end up, in two dimensions, to have only two 
independent components T = Tzz = - T^'^ + 2iT^'^) and f = Tjj = \{T^^ - T^^ + 

2iT^'^), for which the continuity equations read 

d,T{z, z) = dzT{z,z) = 

expressing that T and T are functions only of z and z respectively: T = T{z) and f = 
T{z). The theory separates into two "non-interacting" parts with dependence on z and z 
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respectively. In what follows we shall concentrate on the z part, although all statements 
can be repeated for the z part too. 

The stress-energy tensor is the generator of conformal transformations, in the sense 
that if we take its Laurent expansion 

r(.)=$:A_ ,e. L„ = ^i|T(.).»« (3) 

then the mode L„ generates the infinitesimal conformal transformation 2; — > 2 + ez^^^ 
(analogous expansion holds for Tiz) with modes In particular L_i generates the 

translations in the z direction, Lq + generates the dilations, Lq — Lq the rotations. One 
can write down conformal Ward identities, that, once integrated, give the operator product 
expansion (OPE) for any field with the stress-energy tensor. In particular the OPE of T 
with itself reads as 

r^, ^ c/2 2T{w) dT(w) 

T{z)T{w) = - — - — r- + 7 + — + regular terms 4 

[z — w)^ [z — wy z — w 

The number c is called conformal anomaly or often central charge and plays a central role 
in the following. Eq.(|) implies for the modes (|) the Virasoro (V) algebra 

[L^, Ln] = (m - n)L^+n + — (m^ - m)5m+n,o (5) 

where the constant c appears in the central extension term. The algebra of the full con- 
formal group (do not forget the z part) is then V ® V. We shall use the so called radial 
quantization^ where the time is taken as log \ z\. In this description Lq + Lq is the Hamilto- 
nian, and its eigenvalues must be bounded by below. All the states of a CFT must lie in 
some irreducible representation of the algebra V ® V. The representations we need must 
have Lo and Zq eigenvalues bounded by below, i.e. they must contain a so called highest 
weight state (HWS) |A) for which 

Lo|A) = A|A) , L„|A) = , n>0 (6) 

These representations are known as highest weight representations (HWR). The irreducible 
representations of V are labelled by two numbers; for HWRs these are c and A. We shall 
denote the HWRs of V by Vc(A). For a given theory, c is fixed by eq.(|), then the Hilbert 
space Ti of the theory is built up of all possible representations Vc(A) at fixed c, each one 
with a certain multiplicity: 

H = 0Ar^,A^e(A)®V,(A) (7) 

If a certain Vc(A) ® Vc(A) does not appear, then simply A/a, a = 0. The numbers A/a, a 
count the multiplicity of each representation in 7Y, this implies they must always be non 
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negative integers. They are not fixed by conformal invariance. Constrains on them arise 
from other physical requirements such as locahty ||12|| or modular invariance f^. 

From eq.(^ one easily computes the two point correlator of the stress energy tensor, 
and comparing with eq. (H) it turns out that c is exactly the value of the c{g) function of 
the c-theorem at this particular fixed point. 

Any state \s) = |A,A) in the theory can be put in 1 to 1 correspondence with 
a field through the formula \s) = As{0,0)\0) , where the vacuum |0) is projective (i.e. 
Lq, Lq, L±i, L±i) invariant. In particular the HWS (|D correspond to some fields 0a, a(^! ^) 
that transform under the conformal group as 

0a,a(^,^-)=(^)^(§)%a(^',^-') (8) 
They are called primary fields. Their OPE with the stress-energy tensor is given by 

I {z)(pA ^{w,w) = — — 1 h regular terms (9) 

[z — wY z — w 

Applying this formula to the vacuum and using @ we go back to (P). A basis for all 
the other states of the theory (called secondaries) can be obtained by applying strings of 
Ln,n < to I A). The commutation relations imply 

LoL:;|A) = (A + nA;)Lj:|A) (10) 

Therefore Lq eigenvalues organize the space Vc(A) (often called a module) so that the 
states lie on a "stair" whose 7V-th step (called the A^-th level) has Lq = A + A^ 

states level Lq 



L_3|A),L„2^„i|A),L3_JA) 
L„2|A),L2_i|A) 
L_i|A) 
|A) 

All the fields corresponding to the HWR Vc(A) are said to be in the conformal family [(f) a] 
generated by the primary field (j)A- 



3 
2 
1 




A + 3 
A + 2 
A + 1 
A 



(11) 



3.2 Unitarity 

We said that for a given theory, c is fixed by eq.(||). Conversely, we can ask how many 
theories are there at a certain value of c. The answer to this question for all c's is what is 
called classification of CFTs. 

As a first step towards this classification, one can ask to describe how many HWR 
Vc(A) are there at a fixed c. In general any value of c and A could work, but if we require 



unitarity, i.e. absence of negative norm states, then we can state the following results 11 
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1. for c < no representation is unitary 

2. for < c < 1 the following set of Vc(A) is unitary 

6 

'=1-^^ . P=2.3.4.... (12) 

A^A^^ KP+'ipPf -1 l<,<r<p-l r,seZ (13) 
4p(p + 1) 

3. for c > 1 all representations are unitary 

Therefore a theory containing negative A's or c is automatically non-unitary. Unitarity 
is an essential requirement in string theory. Also many statistical systems enjoy it, but 
there are well known cases (percolation, Lee- Yang edge singularity) where unitarity does 
not hold (i.e. the Hamiltonian is not real). 

3.3 The OPE-algebra and correlation functions 

All fields in a CFT are expected to obey a closed algebra under OPE 

A{z,z)A,{0,0) = 5^4(2,^)Afc(0,0) (14) 

k 

where i,j, k run over all fields (primaries and secondaries) of the theory. Conformal invari- 
ance (via eq.(|^)) reduces the problem to the knowledge of the OPE-algebra among primary 
fields and puts constraints on the form of the functions B^j {z, z) . The OPE-algebra of pri- 
mary fields reads as 

Uz,z)4>,{0,0) = 5^ [0,(0,0)] (15) 

k 

where now the indices run over all primaries of the theory and [0,(0,0)] means 

contribution from the whole conformal family [0,], which can be seen as an expansion 
over all secondaries of [0,], whose coefficients are also (in principle) fixed by conformal 
invariance. The only objects that remain unfixed are the structure constants C^j. Were 
these known, one could reduce via iterative applications of OPEs ([15|) all the correlators 



among primaries to 2 and 3 point functions, which are fixed by projective invariance 

(0i(zi, ^1)02(2:2, Z2)) = 5i2zl2^zlf' (16) 

(01(Z1, -Zr)Uz2. ^2)03(^^3, ^3)) = Cl^zJ^zJfzlfzlfzlfzlf (17) 

where Zah = Za — Zb, 'jab = Ac — Aa — Afo with a ^ b ^ c and a,b,c = 1, 2, 3. Correlators 
among secondaries could be reduced to correlators among primaries via the OPE (||). So, 
at least in principle, all correlators in the theory are computable: the theory is exactly 
solvable. Unfortunately in most cases we do not know the constants C^j. Constraints on 
C^j come from the requirement of associativity of the OPE-algebra, which is equivalent to 
ask duality of the 4-point functions. So in general C^j can be computed if we know the 
4-point functions. 
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3.4 Null vectors and minimal models 

It can happen that in a certain HWR a secondary state |x) (to be called null vector) 
satisfies Lo\x) = (A + N)\x) and Ln\x) = 0, n > 0. This means that |x) behaves like 
a primary state. This apparent contradiction is solved by saying that the HWR is not 
irreducible. Taking \x) = (that can always be done consistently as it is easy to prove 
that {s\x) = for all \s) G HWR, including |x) itself) we get rid of the representation 
embedded in the HWR. If there are more than one null vector, we repeat this procedure 
until we have "cleaned" the representation from all null vectors. What remains is a true 
irreducible module Vc(A) that can be used to build up a CFT. 

When null vectors appear, they give constraints in the form of partial differential equa- 
tions on the correlators of primary fields. Indeed, consider a correlator {x{z)4>i{zi)...(f)n{zn)) 
where (j)i{zi) are primary fields and xi^) is a null vector in the conformal family [(f)]. 
Then x(z) is a secondary of (f){z) and eq.(|) (or more precisely the conformal Ward iden- 
tity obtained by it) implies the existence of a suitable differential operator V such that 
{x{z)(f)i{zi)...(f)n{zn)) = V{(j){z)(f)i{zi)...(f)n{zn)) ■ On the other hand x = 0, so that we have 

P(0(2;)0i(2i)...0„(2„,)) = (18) 

which is a differential equation to be satisfied by the correlators among primary fields. This 
constraint can be so powerful to select a finite number of primary fields under which the 
OPE-algebra closes. In this case only a finite set of Vc(A) are used to build up the model. 
We shall speak in this case of a minimal model of V. Minimal models exist for c < 1 only 
(there is a theorem by Cardy 0| that prevents from the possibility to construct a model 
with finite number of Vc(A) for c > 1). More precisely unitary minimal models can exist 
only for the values ofc = |,^,|,|,... given by the formula ( |l2l) and they can be built up 
using only the Vc(A) representations such that A is contained in the Kac-table given by 
eq-dl). 

The differential equations ( p!^) reduce, in the case of 4-point functions and after mapping 
zi — > 0, Z2 — >■ 1, -^4 — ^ oo, to ordinary fuchsian differential equations in z = Z3. In the case 
of minimal models a solution to these equations can be given [l^. This means that the 

can be computed, and so exact solvability of the minimal models is ensured. 



3.5 CFT on a cylindrical geometry and FSE 



The mapping of the complex plane into a strip of thickness L with periodic boundary 
conditions is a conformal transformation, realized by the mapping 



22L, 

e L 



e L 



(19) 



with u = T + ia and u = t — ia, where r g] — 00, +oo[ is the cylinder "time" and a G [0, L[ 
the cylinder "space". Under conformal mappings, the stress energy tensor transforms as 



u 



Dm 



dzU 
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Consequently, the modes of T on the plane map on those on the cylinder according to 



inu 

7 fc- 



i.e. they become Fourier modes on the cylinder. The radial quantization in z briefly 
mentioned above becomes in this cylindrical context a usual time ordering quantization. 
Periodic boundary conditions mean T{u + L) = T{u). The physical interpretation of modes 
on the cylinder is that total energy and momentum of the system are given by the zero 
modes 

E = -{Lo + Lo-^) , P = -(Lo-Lo) 

Therefore a state |A, A) on the cylinder represents a state of given energy and momentum 

27r - c 271 - 

In particular the vacuum of a unitary theory (A = A = 0) has energy 

E =- — 

-'-'vac „ T 

bL 

given by the Casimir effect of being on a finite and periodic geometry. We see that the 
constant c has the physical meaning of measuring this Casimir effect. 



3.6 Integrability of CFT 

CFT is an integrable QFT. This can be shown by explicitly constructing an infinite set of 
conserved charges. Define the charges 

Is+i = j daTs{u) = ^ dzTs{z) 

where the currents Tg{u) are time ordered (in r time) polynomials in T and its derivatives 
at level s in the Verma module of the identity operator. One can show that if s is odd, all 
the Ts fields at level s can be arranged to be total derivatives of some Ts-i and therefore 
their integrals Ig+i are trivially zero. Nontrivial integrals of motion generated by T exist 
only for even s. At a given level there can be more than one monomial in T and its 
derivatives, the most general current will be a linear combination of these. If we require 
that the corresponding charges commute with all the ones at smaller values of s then the 
coefficients turn out to be fixed in a unique way, thus giving a sequence of currents, one 
for each even level s, that are in involution, i.e mutually commute. This actually defines 
an integrable QFT. Here are the first few currents of this kind. 

T2{u)=T{u) 
Uu) =: Tiuf : 
niu) =: T(n)3 : +^ : (9„T(n))2 : 
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etc... Of course a similar sequence of integrals Is+i holds also for T. Note that the first 
current is T itself. Therefore, energy and momentum are conserved quantities belonging 
to this set, and build up a vector of spin 1. All the other charges, as they commute with 
the Hamiltonian, are also conserved charges in an infinite sequence of higher and higher 
odd spins. 

Of course, if the model possesses higher symmetries (like in ly-algebras or in WZW 
models) there can be other sets of conserved currents and the sequence of spin of the 
conserved charges could be different. 

3.7 c=l Conformal Field Theory 

A particular example of CFT that will be used frequently in the sequel of this review is 
the case of the so called massless free boson compactified on a target space of a circle of 
radius R. Here we give a brief summary of such a c = 1 CFT. The Lagrangian is taken to 
be 



C=^ I d^ip{a, r)9^y.(a, T)da , a G [0, L] , 



(20) 



where L is the spatial volume (i.e. the theory is defined on a cylindrical space-time with 

circumference L). The superselection sectors are classified by the U{1)^ x U{1)^ Kac- 
Moody symmetry algebra, generated by the currents 

J{z) = idzf , J{z) = idzf ■ 

where we have mapped the cylinder back to the plane according to eq.(p!9|). The left/right 
moving energy-momentum tensor is given by 



oo ^ oo 

T{z) = —dzipdzip = ^ LkZ"'-^ , f{z) = —dz^dzif = 



fc=— oo A;=— oo 



The coefficients L„ and L„ of the Laurent expansion of these fields generate two mutually 
commuting Virasoro algebras. If the quasi-periodic boundary conditions are required for 
the boson 



^[a + L,t) = (p{(T, t) + 2TTmR , m G Z , 

then the sectors are labelled by a pair of numbers (n, m), where ^ {n is half integer because 
of the locality, see later) is the eigenvalue of the total field momentum ttq 
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and m is the winding number, i.e. the eigenvalue of the topological charge Q defined by 



1 '•^ 



In the sector with quantum numbers (n, m), the scalar field is expanded in modes as 
follows: 



y9(a,r)= 0(2;) + 0(2), 

0(2)= \^o-ip+^ogz + iY,k^Q\akZ- 



■k 



- |<^o - log 2; + z iafc2; , 

where the left and right moving field momenta p± (which are in fact the two U{1) Kac- 
Moody charges) are given by 

P± = ;^±2^^- (21) 
The Virasoro generators take the form 



k=—oo k=—oo 

where the colons denote the usual normal ordering, according to which the oscillator with 
the larger index is put to the right. 

The ground states of the different sectors (n, m) are created from the vacuum by the 
(Kac-Moody) primary fields, which are vertex operators of the form 



V(n,m){z,z) =: expi{p+(j){z) + p^(j){z)) : . (22) 

The left and right conformal weights of the field V(„,m) (i-e. the eigenvalues of Lq and Lq) 
are given by the formulae 

A± = ^. (23) 

The Hilbert space of the theory is given by the direct sum of the Fock modules built over 
the states 



|n, m) = V(„,„)(0,0) \vac) , (24) 
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with the help of the creation operators a_fc , a_fc A; > 0: 

H = ^ {a-fci . . .a-kj,a_i^ . . . a_/Jn, m), fci, . . . kp, k, . . . Iq e Z+} 

(n, m) 

The boson Hamiltonian on the cylinder is expressed in terms of the Virasoro operators as 

Hcft = y{^o + Lo-^) , (25) 
where the central charge is c = 1. The generator of spatial translations is given by 

P = ^ (Lo - lo) . (26) 
The operator Lq — Lq is the conformal spin which has eigenvalue nm on the primary field 

^n, m) • 

One can also introduce twisted sectors using the operator T that performs spatial 
translations by L: x —>■ x + L. The primary fields V(n,m) as defined above satisfy the 
periodicity condition TV(^n,m) = ^(n,m)- If the more general twisted boundary condition 
labelled by a real parameter u is required 

'T'V^n,m) = exp (il^Q) V^n,m) , 

then it is possible to generate superselection sectors for which n E Z + 

It is important to stress that a particular c = 1 CFT is specified by giving the spectrum 
of the quantum numbers (ra, m) (and the compactification radius R) such that the corre- 
sponding set of vertex operators (and their descendants) forms a closed and local operator 
algebra. The locality requirement is equivalent to the fact that the operator product ex- 
pansions of any two such local operators is single valued in the complex plane of z. This 
condition, which is weaker than the modular invariance of the CFT, is the adequate one 
since the theory is considered on a space-time cylinder and we do not wish to define it on 
higher genus surfaces. 

By this requirement of locality, it was proved in fl^ that there are only two maximal 
local subalgebras of vertex operators: Ab generated by the vertex operators 

{V{n,m) ■ n, m e Z} , 

and Af generated by 

{V(n,m) ■ n e Z, m G 2Z orn G Z + ^, m G 2Z + 1} . 

Other sets of vertex operators can be built, but the product of two of them gives a nonlocal 
expression. 
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4 Perturbed Conformal Field Theory 



In the previous sections we have seen that a QFT can be thought as a RG flow in the space 
of actions, and its UV and IR points correspond to scale invariant CFT. We have then 
concentrated on CFT and explored its properties. Armed with this new insight we may 
wonder if the CFT data are able to predict some piece of information about the theory out 
of the fixed points. 

We may think to define a QFT as a deformation of a critical CFT by some operators 
[T^ , i.e. to assume as action the following expression 

n „ 

S = ScFT + Y,^i / (fx^i{x) (27) 
i=l 

Of course, the class of QFT2 is larger than the one described by this sort of action. Nev- 
ertheless, this class of the so called Perturbed Conformal Field Theories (PCFT) is very 
important because most of the predictions can be done working in this realm and many 
applications to the vicinity of critical points in the theory of critical phenomena can be 
described by this class of actions. 

Lorentz invariance in 1+1 dimensions is equivalent to rotational invariance in the 2 
dimensional space once a Wick rotation is performed. Hence, to achieve Lorentz invariance 
we must require that the fields $i are scalars under 2 dimensional rotations, i.e. their 
conformal spin Aj — Aj must be null, which implies Aj = Aj. Their anomalous dimensions 
are therefore 2A.j. 

4.1 Relevant, Irrelevant and Marginal Perturbations 

Let us consider the simple case of one single perturbing field $, of left conformal dimension 
A. The coupling A is therefore dimensionful, with scaling dimension y = 2 — 2A. Renor- 
malization is needed [l^, and can be achieved by introducing a dimensionless coupling g, 
a mass scale /i and a renormalized field ^{x,g) by the relations 

g = l2~yZgX , ^x) = y^^x,g) 

where Zg and Z$ are renormalization prefactors collecting all infinite contributions. 

(7 = is a fixed point, because the interaction perturbing term disappears and one is 
left with the pure Soft- The field $ = || must have a limit in the unperturbed CFT, say 
$0, with 

r(0)$o = 7(0)$o = 2A<l>o (28) 

We are asking which are the conditions on the field $ that guarantee that the fixed point 
at (7 = is an UV one, from which the RG flow actually gets out and, increasing t goes 
to some IR destiny. To have this, g must increase towards an IR flxed point g*, and the 



17 



anomalous dimension of $ changes according to the already seen formula 



r(^)$=(2-g)$ (29) 



Putting together eqs.(p^) and ( |29D we can compute the first order of the beta function: 
Pig) = yg + 0{g'^). The nature of the fixed point (UV or IR) is chosen by the sign of the 
first derivative of f3{g), as we saw in section 2. y = /5'(0) > only if A < 1. 
Fields are classified with respect to the RG group as 

• relevant if A < 1. If such a field perturbs a conformal action, it creates exactly the 
situation described above, i.e. the theory starts to flow along a RG trajectory going 
to some IR destiny. 

• irrelevant if A > 1. Such fields correspond to non renormalizable perturbations 
which describe the neighborhood of IR non trivial fixed points. Usually one speaks 
of an attraction field better than a perturbing one. We shall not deal with this case 
in the following, but the interested reader may consult, for example [|l^] to see some 
possible applications of this situation. 

• marginal if A = 1. These are further distinguished by the behaviour of higher 
derivatives of the beta function 



— if they stay marginal at all orders in perturbation theory, they are called truly 
marginal, and describe deformation of a CFT leading to another CFT with the 
same value of the central charge. They identify a submanifold of CFT in the 
space of actions. 

— if instead the anomalous dimension becomes smaller than 2 at higher orders in 
perturbation theory, the fields are called marginally relevant. They provide a 
perturbation of CFT defining a non scale invariant theory, but the perturba- 
tion theory is much more complicated and in particular the smoothness of the 
deformation of the Hilbert Space is no more valid. This sort of UV limit, that 
Liischer defined as singular many years ago, is typical of the asymptotically free 
theories and hence very interesting for the properties that such theories share 
with QCD. 

— finally there can be some marginally irrelevant fields, defining a non renormal- 
izable attraction to an IR no trivial point of particular theories 



4.2 Conformal Perturbation Theory 

In a PCFT it is natural to define a Conformal Perturbation Theory (CPT). This is a 
perturbation theory where the unperturbed part of the action is not taken, as usual, to 
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be the free part, but better the CFT at UV. A general correlator {X) of a string of fields 
X = Ai...An expresses as 



{X) = J P(^Xe-^[^l = J p<^Xe-^^^^-^-^'^'^*(^) 

and expanding in powers of A one can express {X) as a series of conformal correlators (in 
principle computable by CFT techniques) with the insertion of the perturbing field 

(X) = f; f {X^xi)...<l>{xk))cFTdW..d'xk 

The possible IR divergences of such an expansion are cured if we put the theory on a finite 
geometry, e.g. on a periodic cylinder as above. A full treatment and detailed expressions 
for the coefficients of the CPT on the cylinder are given in CPT turns out to be UV 
finite if y > 1 i.e. A < |. For < y < 1 only a finite number of terms are divergent: 



the action (P7|) is super-renormalizable in this case. Once the divergences cured (usually 
analytic continuation in y does the job cleanly) the series are convergent within a finite 
radius (which is not the case in general for perturbation series of higher dimensional QFT). 
It can even happen that the IR point g* falls into the convergence region, like it has been 
shown in |l^ for the minimal models perturbed by $1,3 at large p. In this case, one can 
do predictions on the RG fiow and follow scaling functions of physical quantities by use of 
the CPT. In general, the radius of convergence is too small to contain the IR region of the 
RG flow. Nonperturbative effects take place and one has to resort to alternative methods 
of calculation to study scaling functions in the IR regime in order to relate them to the 
particle content of the theory. 

For y = one exactly gets the case of an asymptotically free QFT. The perturbing oper- 
ator is marginally relevant and the theory is renormalizable but not super-renormalizable, 
i.e. there are in general infinities at each order, but they are cured by the insertion of a 
finite number of counterterms. Finally, for y < the theory is non- renormalizable and 
in general there are new infinities and new counterterms to be added at any order, thus 
loosing predictability. However, in the case of integrable theories, one can go a bit further 
T6|| in restricting the form of counterterms and getting some physical information on the 



attracting field (these theories better describe the attraction to an IR point of a RG fiow 
by some irrelevant operator). 



4.3 Integrability in PCFT 

We have seen that any CFT is an integrable theory and possesses an infinite set of mutually 
commuting charges. The perturbation ( |27D usually destroys such conservation laws. It can 
happen, however, that a subset of the local charges of section ^]6| keeps conserved because 
the modification of the current due to the perturbation can still be arranged as total 
derivative, thus recovering the continuity equation even off-criticality. In more precise 
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words, the equation dzTs{z) = valid in the unperturbed CFT is no longer valid off- 
criticality. It is deformed by contributions coming from the perturbation series 

d,T, = XRi'l, + X'Rfl, + ... + A^/2S (30) 

The fields -R^J*! must have, in the super-renormalizable limit in the unperturbed 

CFT, defined by some field in the algebra of fields of the CFT, with conformal dimensions 
{A^"'\ A^")) appearing in the Kac table of primaries or secondaries. Dimensional balancing 
of eq.(l30|) gives 

A(") = s - n(l - A) A(") = 1 - n(l - A) (31) 



where A is the dimension of the perturbing field $. This implies that the series ( pOD 
actually truncates at some integer as conformal dimensions in a CFT are bounded by 



below. In general there is only one field in the Kac table compatible with (|3lD appearing in 
(^). In this case we speak of a no-resonance perturbation. Concentrating on this simple 
case, where only R^^li appears with conformal dimensions (s — 1 + A, A), we see that such 
field must be a secondary of the perturbing field $ of spin s — 1. It is obtained by applying 
only left modes L_„ to If the combination of modes can be put in the form of 

L_i ■ (any combination of strings of Iv_„'s) 

then, as L^i = dz, one can see that XR[^1i = dz'ds-2 for some field Qs-2 also in the 
secondaries of $ but one level below. Then eq.(|30|) boils down to a genuine continuity 
equation 

d-zTs = dzQs-2 

guaranteeing a conservation of the charge 

Is+i = j) dzTs{z, ^) + y dzQs-2iz,z) 
If such a situation is realizable or not depends on the nature of the perturbation $. 



A. Zamolodchikov [Tjl has given a sufficient but not necessary condition to ensure the 



presence of a conserved current, based on the analysis of the dimensions of subspaces of 



the Verma modules of $ and the identity. For details see the original article ||T4| . 



Integrability of a two-dimensional massive QFT has strong implications on the proper- 



ties of its S-matrix 32 



1. the number of particles is conserved in the scattering process, there is no particle pro- 
duction. Only exchanges of internal quantum numbers are allowed among particles 
in a multiplet; 

2. the final set of momenta coincides with the initial one: {pi, ....,Pn}m = {Pi, ■■■■,Pn}out'-, 
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3. the S-matrix for a process with n incoming and n outgoing particles factorizes into S- 
matrices of 2 incoming and 2 outgoing particles. To guarantee that the order of this 
factorization is irrelevant, the 2 particle S-matrices have to satisfy a factorization 
equation known as Yang-Baxter equation (here 1,2,3 label the particles and % = 
9i — 9j, 9i being the rapidity of particle i) 

S^!^ (012) 3^:^(013)3^:^(623) = ^^i|(^12)^S(^13)^£^(^23) 

The factorizable S-matrix of an integrable PCFT can often be conjectured by the knowledge 
of integrals of motion, that constrain the bootstrap equations and by the Yang-Baxter 
Equation itself, that imposes severe constraints on the matricial form of 3^f{9). However, 
this conjecture has to be verified against independent checks. The particle description of 
the Hilbert space must be linked with the conformal description valid in the vicinity of the 
critical point, otherwise one cannot assert that the two theories defined as perturbation of 
a CFT and as factorized scattering theory are the same. 

4.4 Scaling functions 

One possible bridge between the two formulations - perturbed CFT versus Factorizable 
Scattering Theory - is given by FSE, more specifically by the use of the so called scaling 
functions. We have seen that the energy of a state \i) on a cylinder of circumference L in 
CFT is given by 

Ei = -^ with Ci = c-12(A, + Ai) (32) 

It is obvious that when we move outside the critical point such energy level has a different 
dependence on L. However, the overall 1/L dependence is fixed by dimensional reasons, so 
all the variation of the dependence must be confined in a dimensionless function that we 



are free to normalize in such a way that in the limit L — * it reproduces eq. (|32D . Being 
dimensionless, such a function can depend on L only in a dimensionless way, i.e. we are 
forced to introduce a mass parameter m to compensate for the physical dimension of L. 
The dependence will be on a dimensionless parameter / = mL. This is very reminiscent 
of what already illustrated in the introduction and the discussions made there about the 
UV and IR limits recovered for / — > and / — > oo respectively of course apply here as well. 
They show that for each state \i) in the theory there is a function Cj(Z) attached to this 
state in such a way that its energy is 



Ei(L) 

such that 



T^Ci{l) 

6L 



lim;^oCj(0 = Ci = c — 12(Aj + Aj) i.e. it reconstruct conformal data in the limit 
/ ^ 0; 
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• it can help, in a manner that will become apparent later, to extract scattering data 
in the limit / — oo 

The forced introduction of a mass scale parameter should not surprise: we are actually 
going out of a critical CFT invariant point into a theory which is not scale invariant. It is 
therefore natural that a mass scale should be introduced. 

To such functions Ci{l) the name of scaling functions is usually given. The problem 
to bridge between the PCFT and the Factorized Scattering formulations of an integrable 
QFT2 is recast into the computation of such scaling functions. 

Of course, one way to compute them, also in non-integrable theories, is to resort to CPT. 
However, as already commented, the validity of CPT is normally confined to a small radius 
of convergence that does not allow to make contact with the IR region where scattering 
theory is apparent. This justifies the seek for the development of nonperturbative methods. 

• The Truncated Conformal Space Approach (TCSA) consists in diagonalizing the 
truncated {E < Ecut) Hamiltonian Hq'ft + Vpert numerically. It is non-perturbative, 
and applicable even to non-integrable QFT. However, it is affected by truncation 
errors that increase with / and, being totally numeric, does not allow for any analytic 
control of the scaling functions. 

• The Thermodynamic Bethe Ansatz (TBA) [3^ implements the thermodynamics of 
a gas of particles governed by the given S-matrix. The free energy can be exactly 
computed in terms of a set of coupled non-linear integral equations. Exchanging 
the role of time and space, such free energy can be reinterpreted as a FSE Casimir 
vacuum energy where the role of L is played by the temperature. Recently also 
excited states have been accessed with this method [3^, that has been very popular 
for the last ten years. The treatment of the vast literature on this method is out of 
the scope of this review. 

• In the next sections we shall illustrate an alternative method to obtain scaling func- 
tions for an integrable model starting directly from its definition as a QFT regularized 
on a lattice. This method is know in QFT as Destri De Vega Non-Linear Integral 
Equation (NLIE) from light cone Bethe Ansatz ||2^, but similar equations have also 
been derived in Condensed Matter Theory ||35| . 

4.5 Sine-Gordon / massive Thirring models as PCFT 

The sine-Gordon model, which has been well known at the classical level for the late fifty 
years, plays also an important role in quantum theory, thanks to its particular properties 
of non-linearity and integrability. It has been successfully applied to very different sectors 
of Mathematics and Physics, ranging from partial differential equation theory to parti- 
cle physics or solid state physics. Recent applications of the classical model are related 
to nonlinear optics (resonant dielectric media) and optical fibers, magnetic properties of 
polymers, propagation of waves in crystals, etc... Interesting applications of the quantum 
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model are related to Kondo effect and to the thermodynamics of some chemical compound, 
as Copper-Benzoate Cu(C6D5COO)2-3D02 ||l9l. At the same time, the quantum theory 
shows a phenomenology that is similar to the Skyrme model used before QCD era to 
describe barions and strong interactions. 

The most relevant properties of the model are 

• at a classical level, all the solutions of the equations of motion are known (exact 
integrability via inverse scattering method) 

• the classical solutions describe solitons, antisolitons and bound states (breathers); in 
a scattering process this solutions are transparent (it is the mathematical meaning 
of "soliton") 

• it admits, both at a classical and at the quantum level, a countable infinite set of 

conserved charges 

• the quantization of the theory describes an interacting particle with its antiparticle 
and, in a certain (attractive) regime, bound states 

• the S matrix has been exactly determined; only elastic scattering processes can take 
place (i.e. no particle production), that is the quantum analog of the classical trans- 
parency of solitons 

The minkowskian Lagrangian of sine-Gordon theory is given by 



where $ denotes a real scalar field, while that of the massive Thirring theory is of the 
following form: 



C^Th = j {^^{i^.d^ + m^)^ -U>Y'^^lu'^)dx, (34) 

describing a current-current selfinteraction of a Dirac fermion It is known that the two 
theories are deeply related provided their coupling constants satisfy 
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47r 1 + gj-K 

For comparison with the Destri-de Vega nonlinear integral equation, it is important to deal 
with cylindrical sine-Gordon and massive Thirring, i.e. the integrals in (|33| , |34|) must be 
taken in the interval [0, L\. 



The cos term in ([33|) can be considered as a perturbation of the c = 1 free boson 
compactified on a cylinder, as described in section p77| . Similarly the massless (mo = 0) 
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Thirring model is a c = 1 conformal field theory and the mass term plays the role of a 
perturbation. From Coleman's paper pO|] it is known that correlation functions of the 
perturbing fields and : cos/3$ : are identical, then both models can be considered as 
the perturbations of a c = 1 compactified boson by a potential V : 



SsG/mTh = ScFT + V , V = \j {V(i^o){z, z) + V(^^i^o){z, z)) (f 



X 



(35) 



which is related to the bosonic Lagrangian (|33D by the following redefinitions of the field 
and the parameters: 



— A = 



^ = VA7r^, R = ^, >^ = WM- (36) 



For later convenience, a new parameter p can be defined by 

/^^ 1 

P = ^^^^ = ■ (2^) 

The point p = 1 (i.e. g = 0) is the free fermion point, corresponding to a massive Dirac free 
fermion. The particle spectrum of Sine-Gordon for p > 1 is composed by the soliton (s) 
and its antiparticle, the antisoliton (s). It is known as repulsive regime because no bound 
states can take place, p < 1 is the attractive regime, because s and s can form bound 

states that are known as breathers. The values p = t y k = 1,2, . . . are the thresholds 

k 

where a new bound state appears. The potential term becomes marginal when = Svr 
which corresponds to p = 00. The perturbation conserves the topological charge Q, which 
can be identified with the usual topological charge of the Sine-Gordon theory and with the 
fermion number of the mTh model. 



Mandelstam ||2T| showed that a fermion operator satisfying the massive Thirring equa- 
tion of motion can be constructed as a nonlocal functional of a pseudoscalar field (boson) 
satisfying the sine-Gordon equation. But the fermion and the boson are not relatively local 
and then do not create the same particle (the two theories are not equivalent). 

The difference between them is that they correspond to the perturbation by the same 
operator of the two different local c=l CFTs Ab and Aj as in |3]^. The short distance 
behaviour of the sG theory is described by the local operator algebra At, while the primary 
fields of the UV limit of mTh theory are Af. 

Note that the two algebras share a common subspace with even values of the topological 
charge, generated by {V(„ m) : n G Z, m G 2Z}, where the massive theories described by 
the Lagrangians (l3^ ) and (|34|) are identical. The well known proof by Coleman about the 
equivalence of the two theories ||2^ holds exactly in this subspace. Figure |^ shows the four 
sectors where all the vertex operators live. 
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Figure 3: The family of vertex operators V(n,m) with n £ Z/2 and m £ Z. Sector I is the common 
subspace. I and III are Aj, that defines the UV behaviour of massive Thirring; I and II are Ab 
that is the UV of sine-Gordon, IV is a sector of non mutually local vertex operators. 

5 Light-cone lattice regularization of Sine-Gordon the- 
ory 

In this section v^^e present a lattice regularization of the Sine-Gordon model which is par- 
ticularly suitable to study FSE. It is v^^ell known to lattice theorists that there are many 
actions on a lattice that all discretize the same continuum theory. This means that there 
are many possible regularizations of the same theory. Lattice researchers are used to choose 
the action on the lattice possessing the properties that best fit their calculational needs. In 
the present context what we would like to keep on a lattice discretization of Sine-Gordon 
model is of course the property of integrability. The light-cone lattice construction we are 
going to illustrate is a way (not the unique!) to achieve this goal. 

5.1 Kinematics on light-cone lattice 

It is a usual way to regularize quantum field theories by defining them on a space-like 
"Hamiltonian" lattice (where time is continuous and space discrete) or space and time- 
like "Euclidean" lattice (when both space and time are discrete). In statistical mechanics 
this is not just a regularization method but can be a right microscopic way to describe 
physical systems. In two dimensions, the most known approach is to define a rectangular 
lattice with axis corresponding to space and time directions and associate to each site an 
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interaction depending only on the nearest neighboring sites. In this case the partition 
function can be expressed in terms of a transfer matrix. 

In what follows, a different approach |^] is adopted: Minkowski and Euclidean space- 
time can, in fact, be discretized along light-cone directions. Light-cone coordinates are: 

x± = X ±t 

and the choice 

M. = {x± = -^n±, n± G Z} 
V2 

defines a light-cone lattice of "events" as in figure || (a) . They are spaced by a in the 
space and time directions and by a/V2 in light-cone directions. At every event P G A4 a 
double light-cone (in the past and in the future) is associated and only events within this 
light-cone can be causally connected (see fig. ^ (b)). Then, any rational and not greater 
than 1 speed is permitted for particles, in an infinite lattice. The shortest displacement 
of the particle (one lattice spacing) is realized at speed ±1 and corresponds, from the 
statistical point of view, to nearest neighbors interactions. Smaller speeds can be obtained 
with displacements longer than the fundamental plaquette, and correspond to high order 
neighbors interactions. In quantum field theory, these are nonlocal interactions. In the 




X 



(a) (b) 

Figure 4: (a) light-cone lattice; (b) double light-cone emanating from a point and high-order 
interactions (dashed) 

following, only the local case (nearest neighbors) is treated. The nearest neighbors of the 
event P are the four nearest points in the light-cone directions. This implies that particles 
can have only the speed of light ±1 and are massless. They are called right-movers (R) 
and left-movers (L). 

In this case, it is possible to introduce a useful language for connection with statistical 
mechanics associating a particle to a link. Consider the two links in the future that come 
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out from an event P. Particles R and L in P, by definition, are respectively associated to 
the right-oriented link and to the left-oriented link. In this way, the state of a link is defined 
to be the state of the point where it begins (also the opposite choice, of connecting a link 
with the site where it ends, can be done; it is simply a matter of convention) . For example, 
if in a point O there is a particle R, one tells that the "right-oriented" link outcoming from 
it is occupied by R. This correspondence of points and links is possible because only local 
interactions are assumed, and it is useful because the counting of states is simpler. In the 
following, the lattice is assumed of a finite extent L = aN in space direction {N is the 
number of sites, counted in the space direction), with periodic boundary conditions, but 
infinite in time direction. In this way a cylinder topology is defined for space time. The 
Hilbert space of states in an event P is the tensor product 

of R and L space of states. The fact that particles can be classified in left and right does 
not mean, in general, that the two dynamics are independent, as it happens in CFT. 

Call |aLi,affi) the generic vector of a basis of Hi where i = 1,...,N labels the sites. 
The notation 

|a2i-i,a2i) = |aLi,affi) 

is useful and not ambiguous (even number refers to right, odd number refers to left). The 
total Hilbert space is: 

N 

Hn = Hi 

1=1 

and a basic vector can be represented by 

"2) ® ••• ® \a2N-l, C12n) = «2, «2Ar) G Hn- 

Note that in a sites lattice, due to light-cone, 2N labels are required. The state 
|ai, 02, a2Ar, t) is better thought to lie on links separating two lines of events. The 
time slices t, t + |, t + a are taken as represented in fig.||. From this figure it is clear that 
two different types of evolution operators can be defined, depending on the initial state: 

f/+ |ai, 02, a2N, t) = \a[, a'2, a^^, t + a/2) 

(38) 

U- \a[,a'2,...,a'2N,t + a/2) = \a'l,a'^, ...,a'-^^,t + a) 

where the initial states are chosen as in figure ^ They can be pictorially represented by 
the links that have to be added to the time slice t to get t + f (for U+) and to t + | to get 
t + a (for U^) 
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Schrodinger form of equations of motion is used, for a state in Hilbert space.Q The 
global time operator can be chosen as 

U = U+U- or U' = U.U+ 

depending on the initial state. For a consistent quantum theory, both these operators must 
be unitary. This is guaranteed if the assumption U\.U+ = Ul_U- = 1 is made, that is the 
elementary operators themselves must be unitary. 

Another operator plays an important role and is defined as follows (the states are at a 
certain fixed time): 

V |q!i, 02, a-ZN, t) = \(y2N, Oil, a2Ar-i, t) (39) 

It corresponds to an half-space shift in the space direction, with exchange of right and left 
states (see the figure ||). 

Two applications of V give a shift by an entire lattice spacing, then V"^ is the lattice space 
evolution operator. Also V is a unitary operator. By using the pictorial representations 
introduced above it is easy to convince oneself that the following relations hold 

[V\ f/±] = 0; U± = VU^Vl (40) 

For example, start by taking f/_. To obtain VU^ one has to apply V after U-, which 
means shift the upper endpoints of f/_ by | towards right. 

^In Schrodinger form, if |a) is a state vector, its time evolution is given by \a,t) = U\a,0) where U 
satisfies motion's equations: U = Te^' '^'**^ (Dyson's series) 



28 



Figure 6: Half-shift operator 



V U_ = 




Now, by applying V"^ to the lower part of this drawing, i.e. shifting the lower endpoints 
by |, one can see that VU^V^ = U^. Analogously, all the other relations can be proven. 

As a consequence of ([40| ) the two principal unitary evolution operators, and U, are 
commuting: 

[V'',U] = [U,V^] =0. 

It is natural to identify them as the exponential of the Hamiltonian operator, and the 
exponential of the linear momentum: 

Y2 ^ g-iaP l4ij 



There are other two important operators, defined as: 

(42) 



Ub. = U+V 



They correspond to one step evolution in light-cone directions. They are commuting and 
give the expressions: 

U = UrUl V = UrUI 

[Ur, Ul] = uIUr = UIUl = 1 

then, using also (^) yields: 

^^ = e-t(^+^), [/^ = e-^t(^-^). (43) 



5.2 Dynamics on light-cone lattice 

A dynamics can be defined by giving all the amplitudes of the different processes that 
can take place on the lattice. The fundamental assumption is that at every site a whole 
process can happen, in the sense that if \aLi(^R)in ^^'^ \(^Li I3r) ^ut ^^e incoming and 
outgoing states in a certain site, they can be considered asymptotic states and the transition 
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amplitude is a (bare) S-matrix element. To each site of the lattice, which is the crossing 
of 2 incoming and 2 outgoing lines, we associate a matrix given by the Boltzmann weights 
of a vertex model 

out {Pr, Pl \aR, aL),^ = i?S>. W (44) 

The parameter A is the so called spectral parameter of the R matrix of vertex models p3[ . 
The R matrix can be pictorially represented by drawing the crossing where it is put and 
labelling the external lines with the 4 indices of the matrix 

Pi P2 

P1P2 

R 

CX J CX 2^ 

The system is then defined via its microscopic amplitudes, formally identified with the 
Boltzmann weights of some lattice model. If the lattice model is integrable, its Boltzmann 
weights, in the form of R matrix, satisfy Yang-Baxter equations, and our bare scattering 
theory also does. This is the way integrability is implemented into the lattice definition of 
the model. 

From the definition (^) of evolution operators one can see that they are realized in 
this context as products of R matrices 

N 

{a,,a2, ...,a2N\ U+ |/?i,/?2, ...,P2n) = H^fclSU^ " (45) 

i=l 

which, in pictorial form, is nothing else than the already seen definition of U^, decorated 
with indices. The operator f/_ can be obtained from (^0|) . The product is on all the 
sites at a given time. The particular choice of the spectral parameter, with a "correction" 
(inhomogeneity) depending on the site will become clear in a moment. 

Apparently, this "phenomenological" approach is quite unusual, because in traditional 
lattice quantum field theory every site is associated with an interaction potential (ex: 0^(i) 
is the potential on the site i), not a whole scattering process. This can appear as a sort 
of "macroscopic" approach, not based on fundamental interactions. But the properties of 
factorizable scattering must be taken into account. Factorization of generic amplitudes in 
2 — > 2 particles amplitudes is a sort of quantum superposition principle and the remarkable 
fact is that between one scattering and the other, the particles are asymptotic ones, that 
means that they are free. Every point contains all the interaction. This is what was 
assumed in the definition of the lattice. Then it is perfectly justified that every site is 
connected with a whole 2^2 particle scattering process.^ 

^In the next paragraphs, it will be shown that, for the particular case of the 6 vertex R matrix, a more 
traditional lattice QFT approach can be formulated, in terms of a fermionic field. 
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5.3 Euclidean transfer matrix 



In order to relate the evolution operators to something computable, we have to consider 
the Euclidean lattice vertex model transfer matrix. Consider a two dimensional Euclidean 
square lattice, with periodic boundary conditions, in both the directions. The links are 
the physical objects of the system. They can be in different states belonging to the vector 
spaces A (horizontal) and V (vertical). In principle, such vector spaces can be different, but 
in the following they will be assumed as isomorphic. At every site we associate a Boltzmann 
weight depending on the four links crossing at this site and on a spectral parameter A (see 
fig. 1^). The simplest case is to take the same Boltzmann weights at all sites, but more 
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Figure 7: (a) periodic (toroidal) square lattice; (b) Boltzmann weight for a site; the labels are: 

a,beA, a,(3eV 



general configurations are possible. For our purposes it will be interesting to consider the 
case with an inhomogeneity Aj at each site, where i is the column index (i.e. all the sites 
on a column have the same inhomogeneity), in the sense that the Boltzmann weight in the 
column i is taken to be w{ab\aP, A — Aj) = Raai'^ ~ -^j)- ^l^o for the boundary conditions 
it is possible to assume more general configurations than the simplest one (i.e. toroidal 
b.c). Assume that between the column and the A^+1 (that is 1) there is a nontrivial 
seam line, in such a way that the Boltzmann weights on the column A^ (with respect to 
the normal ones) are given by e*'^^i?fj^(A — Xn)- This choice is made because only the link 
b crosses the seam line. It implements the so called twisted boundary conditions (of course, 
a; = reproduces the periodic case). 

The object of interest for us is the row-to-row transfer matrix. It is a product of 
concatenated R matrices along a whole horizontal line of the lattice and describes the 
evolution from one time slice to the next in Euclidean lattice. For simplicity, let us insist 
on the pictorial notation 



t{X\{X^) = 



We have omitted the indices on the two rows to avoid useless heaviness of notation. 
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This object is well known in the statistical mechanics literature and for many integrable 
models it can be diagonalized exactly by Bethe Ansatz methods, as we shall illustrate in 
next sections. Unfortunately, it refers to a row-to-row evolution on an Euclidean lattice, 
and not to the case we are interested in of an evolution in tine in our light-cone lattice. Is 
there a way to relate this known object to the operators U± or Ul,r introduced above? 

The trick comes eventually by considering a particular distribution of inhomogeneities. 
Compute the transfer matrix at a value of the spectral parameter, and choose inhomo- 
geneities to be Aj = (— 1)*"''^6. In this case, all the R matrices at odd sites are calculated 
at spectral parameter 29, while those at even sites are calculated at spectral parameter 0. 
The nice fact is that, for all unitary R matrices 



or in graphical form 



Pi P2 

R (0) = 




a 



a. 



The transfer matrix t(6|{(— 1)*+^6}) then turns out to be equivalent to the matrix 
element of the operator U+ taken between two time slices of the light-cone lattice, as it 
appears clearly from this pictorial representation 



t (0|{(-1)'^0}) 



Analogously the operator Ur is realized by t"^(0|{(— 1)*+^©}). Diagonalization of t 
gives therefore the eigenvalues of Ur and Ul which are the exponentials oi E + P and 
E — P respectively. If we are able to compute the eigenvalues of t, we are able to calculate 
the energy levels and the total momenta of states in the theory defined on the light-cone 
lattice. 



5.4 6 vertex model: main results 

The theory developed up to now is general and not referred to a specific model. The 
simplest non trivial case is the 4 x 4 R matrix, corresponding to the choice A = V. As 
shown in ||23|, the most general solution is the so called 8 vertex model. This name means 
that only 8, between the 16 entries of the R matrix, are nonzero. A special case is the 
6 vertex model, for which many results have been obtained in the light-cone description: 
this will be the principal object in the following. The R matrix has the form (lower index 
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are rows and upper index are columns) 



c b 
b c 



(46) 



There is a well known mapping between vertex models and spin chains (see |^), i.e. 
the transfer matrix is the exponential of the quantum Hamiltonian of the chain: 

In the case of 8 vertex model, the Hamiltonian is the XYZ(l/2) chain, while in the special 
case of 6 vertex, is the XXZ(l/2) chain. In what follows, this identification can be useful 
to interpret some facts connected with Bethe Ansatz. The XXZ(l/2) chain Hamiltonian 
is given by: 



N 



and 7 is the anisotropy. The a are Pauli matrices. The total z-component of the spin, 
which is a conserved quantity in this system, will play an important role in Bethe Ansatz. 

It is possible to give a particle interpretation to the same R matrix (^) on the light- 
cone lattice. The simplest case of particles obeying Pauli exclusion principle (fermions) 
and without internal degrees of freedom (color number) is assumed. This means that in 
an event only one particle of type R and one L at most can take place. In other words, 
one link has two states: empty or occupied. At every point there are four links, that 
means 16 possible configurations associated to it. In terms of events, these are the possible 
configurations that connect a point with the nearest neighbors in the future. 

Assume now that only amplitudes that conserve the total number of particles (R+L) are 
non vanishing. This reduces to 6 the permitted configurations, as it is shown in the figure 
||. This is simply the 6 vertex model whose R matrix is written in (|46[). The assumption 









CO, CO. CO, CO. CO^ CO, 

1 2 J 4 3 o 

Figure 8: 6 permitted amplitudes. Dots are particles. 

of integrability for this amplitudes gives the general six-vertex model. The requirement of 
symmetry under parity transformation implies that cua = (^4 and 00^ = 00^. The convention 
adopted in (^61) gives 

b = uj^ = uj4^ c = 005 = Uq. 
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This R matrix can be written in an operatorial form, by defining a lattice chiral fermion 
^A,n, with A = R,L, and n labels the sites. The anticommutation rules are the canonical 
ones: 

{^A,n, ^B,m} = 0, \ i)A,n, m [ = ^AB^nm- (47) 



This fermion has some interesting properties, that are exposed in the paper ||22|, and are 
sketched in the following list: 



1. the R matrix and all the other operators Uany can be written in an operatorial form 
in terms of the fermion; 

2. the lattice Hamiltonian, in the free case ui = uj2 = h = 1, c = 0, can be explicitly 
written; by this, the following dispersion relation can be obtained: E = ±k] this 
is the dispersion relation for a free massless particle; the unusual fact is that it is 
monotonous, then there is no doubling of fermions; this is a consequence of the 
non-locality of the Hamiltonian 

3. the lattice Hamiltonian admits a continuum limit oo and a — > but with 
L = Na fixed; locality is recovered in this limit; the continuum equations of motion 
are those of the massive Thirring model. The space is compactified on a cylinder. 



In the continuum limit, the massive Thirring model emerges as the field theory charac- 
terizing the scaling behavior of the dynamics on the lattice (remember that L is finite; 
the scaling behaviour is understood in terms of this L). We conclude that the theory 
we have constructed is a lattice regularization of the Thirring model suitable to study its 
integrability properties. 



5.5 6 vertex model: Bethe Ansatz 

The assumptions of unitarity and hermitian analyticity will be taken into account, for the 
R matrix, and this requires that the variables in (^61) must have the specific form: 

a = a('i9, 7) = sinh(z7 — -i?), b = b{'t},j) = smh.{}, c = c('i9, 7) = z sin7 (48) 

The transfer matrix defined by this R matrix can be diagonalized with Bethe Ansatz 
method |^. In terms of the spin chain, this means that there are two operators, usually 
indicated by B^'d) and Ci^'d), whose expression is known, and there is a "reference state"^ 
1^2) such that: 

c{d) \n) = 

^The "reference state" here is only a mathematical object. Physically speaking, it corresponds to the 
ferromagnetic state with all the spins up. 
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and 



(49) 



for appropriate values of is an eigenstate of the transfer matrix. The "appropriate 
values" of i}j can be obtained as the solution of a set of M coupled nonlinear equations, 
that are called Bethe Ansatz equations. In general, the transfer matrix contains all the 
conserved charges, in particular the Hamiltonian. Then the Bethe Ansatz eigenstates are 
also eigenstates of the Hamiltonian. The Hilbert space of the theory and the action of 
conserved charges on it are then exactly known. 

All this computations for the 6 vertex model were performed in ||23|, |2j, the final 
results are written here for the eigenvalues t{{},Q,uj) of the inhomogeneous and twisted 
transfer matrix (here '& is the spectral parameter, 6 the inhomogeneity and u the twist) 



M 



""""""" 

sinh ^[i^-^j- ^] 



M 



and the values of are defined by the set of coupled nonlinear equations called Bethe 
Ansatz equations: 



/ sinhl[^, + 9 + f] sinh ^ - e + f ] 

Uinh 1 + e - f ] sinh ^ - e - f ] 



N 



—e 



2iu) 



M 

n 



sinh ^ ["Oj — 'dk + ^^r] 



L sinh li^j -'dk- «7r] 



where 2N is the length of the chain and the number of sites in a row of the light-cone 
lattice. The 'dj are called Bethe roots, and in principle, can take any complex value. But 
there is a periodicity in their imaginary part: 

,2 



vr 

7 



(50) 



then only a strip around the real axis must be taken into account for the Bethe roots: 



X I 



'27' 27 



(51) 



Moreover, only the range < 7 < vr will be examined. 

The whole spectrum of the theory can be obtained using all the Bethe configurations 
having M < N and 'dj 7^ 'dk for every j 7^ k. In general for a state with M roots the third 
component of the spin of the chain is 



S = N-M, 



(52) 
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because every operator -B(#t) counts as —1 spin. 

This XXZ(l/2) chain has 2 states in every site, then 2^^ states in total. Then the 
energy spectrum is upper and lower bounded. Changing the sign of the Hamiltonian gives 
another permitted physical system. This means that there are two possible vacua. The 
first one is the so called ferromagnetic ground state, corresponding to M = that is the 
reference state \Q) . It has spin S = N. The second one is the antiferromagnetic ground 
state, that can be obtained with M = N and all the roots real. It has spin 5* = 0. 

In what follows, only the antiferromagnetic ground state will be considered, because it 
has one important property: in the thermodynamic limit (A^ — > oo) it can be interpreted 
as a Dirac vacuum (a sea of particles created by B) and the excitations on this vacuum 
behave as particles. 

The energy E and momentum P of a state can be read out by the transfer matrix 
eigenvalues. The final form is: 



qiE±P) ^ ^±^. TT sinh^[zf-9±^,] 

Alsinh^[z| + eT^,]' ^ ^ 



M 

6 2 

3 



Other integrals of motion can be obtained in a similar way by transfer matrix. Note that 
the second term in the transfer matrix expression vanishes because 6(0) = 0. 



6 Nonlinear Integral Equation from Bethe Ansatz 

In this section the fundamental nonlinear integral equation driving sG scaling functions is 
derived. In the literature it is known as Destri-de Vega equation. It was obtained first in 



|26| , pSj for the vacuum scaling function. The treatment of excited states was pioneered in 



38|| and refined in p^, B9|, to arrive to the final form, in pSl pO . 



6.1 Counting function 

It is possible to write the Bethe equations (||) in terms of a counting function Zn{i!}). 
First, introduce the function^ 

sinh -^(iniy + •&) 
smh ^^(zTTZ/ — V) 

The oddness on the analyticity strip around the real axis defines a precise choice of the 
logarithmic branch. The counting function is defined by 

M 

Z^i^) = N[<f)y2{^ + e) + 01/2(^9 -Q)]-Y1 - ^k) + 2u (54) 

k=l 



''In the following we use the parameter p = ^ — 1. It will later turn out to be exactly the same as the 
one introduced in (p7|). 
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The logarithm of the Bethe equations boils down to the simple condition 

Z;v(^,) = 27r/,- , eZ + i^, 5= (M)w2 = (iV-^)w2 e {0,1} (55) 



The number Ij plays the role of a quantum number for the Bethe basic vectors (^Qf ) and it 
must be chosen depending on the value of 6. Notice that 6 and u play a similar role, because 
both produce a shift in the quantum numbers Ij {if u is absorbed in the definition of Ij): 
in the first case the shift is exactly vr, in the second case it is a real (possibly complex) 
number. This means that the variable 6 can be absorbed in oo but the most convenient 
choice is to use them both. 

Observe that Bethe roots can be obtained as zeros of the equation: 

1 + (-l)^e*^^(''^) = (56) 
6.2 Classification of Bethe roots and counting equation 



From Bethe Ansatz it is known that a Bethe state (fl9|) is uniquely characterized by the 
set of quantum numbers {Ij}^^^ , M < N that appear in ([551). Notice that M < N 
means 5" > 0. The values of "dj to put in (|49|) can be obtained solving Bethe equations. It 



is also known that only states with 



are allowed. It is a sort of fermionic character for Bethe states |^ 



Bethe roots can either be real or appear in complex conjugate pairs. In the specific case 
(0), there is another possibility, due to periodicity (|50D: if a complex solution has imaginary 
part Ini'd = |(p+ 1) it appears as a single root (in (H) it produces the left hand side real, 
then its complex conjugate is not required). It is called self-conjugate root. Remember now 
that the maximal number of real roots (M = A^) describes the antiferromagnetic ground 
state. From the point of view of the counting function, a more precise classification of 
roots is required: 



real roots; they are real solutions of ( p5| ) that appear in the vector (^Qf ); their number 
is Mr; 



holes; real solutions of (l55|) that do NOT appear in the vector (^9|) ; their number is 

special roots/holes (special objects); they are real roots or holes where the derivative 
Z'j^^-dj) is negative^; their number is Ns; they must be counted both as "normal" and 
as "special" objects; 

close pairs; complex conjugate solutions with imaginary part in the range < 
llrwdl < 7rmin(l,p); their number is Mq; 



^The characteristics of this type of solutions will be clarified in the next sections. 
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• wide roots in pairs; complex conjugate solutions with imaginary part in the range 
7rmin(l,p) < \Im'&\ < vr2±i; 

• self- conjugate roots (wide roots appearing as single); /m-i? = vr^yi; their number is 
Msc. 

The total number of wide roots appearing in pairs or as single is M^y. The following 
notation will be used (sometimes) for later convenience, to indicate the position of the 
solutions: hj for holes, yj for special objects, Cj for close roots, Wj for wide roots. 

Complex roots with imaginary part larger than the self-conjugates are not required 
because of the periodicity of Bethe equations. A graphical representation of the various 
types of solutions is given in figure ^. 



AlmG 



2y 



w. 



X ^2 



W, 



h, Ree 

4 



Figure 9: The different types of roots and holes and their position in the complex plane, n 
denotes 7rmin(l,p). The upper line at ^ = |(p + 1) is the self- conjugate one. 

An important remark must be done: from the definition of Zn (0) it is obvious that 
only for states without complex roots the fundamental strip for 0i.(t9)), that is the largest 
strip around the real axis without singularities, is the fundamental strip for Zj^. In all the 
other cases the analyticity strip for Z^ is narrower, and depends on the imaginary parts 
of the complex roots. 

An important property follows from this classification: the Z^ function is real analytic 
if a; is a real number 

Z^{r) = {ZNm* (57) 

By considering asymptotic values of (pui^) and Zi^{'d) for ?9 ±oo, it is possible to 
obtain an equation relating the numbers of all the various types of roots. We refer the 
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reader interested in the details of the derivation to |2^, |30|. Here we only mention the 
final result, in the form where the continuum limit oo, a — and L = Na finite, is 

already taken 

Nh -2Ns = 2S + Mc + 2 9{p - 1) Mw (58) 

where 9{x) is the step function: 9{x) = for x < and 9{x) = 1 for a; > 0. Recall that 
5* is a nonnegative integer. In the case of = 0, it turns out that Nh is even (Mc is the 
number of close roots, and is even). 

The most important fact is that the number of real roots does not appear in this 
equation . This fact, together to what will be explained in the next paragraph, allows to 
consider the real roots of particles (Dirac vacuum) and all other types of solutions 

(holes, complex) as excitations on this sea. 

6.3 Non linear integral equation 

Let X be a real solution of the Bethe equation. Thanks to Cauchy theorem, an analytic 
function /(x) on an appropriate strip containing the real axis admits the following repre- 
sentation 

where is a anti-clockwise curve encircling x and avoiding other singularities of the 
denominator, i.e. other Bethe solutions (real or complex). In the region where 0i(?9) is 
analytic, we can use (|59D to write 



Mr+Nh Mr+Nh 



Y: m^-x,)= j2 j ^'^^^^-^^ i + (-i)v^.(.) = 

k=l k=l -^^^k ' \ J (go) 

fdfi^,^ (-l)V^ivMzZ;,(/x) 

The sum on the contours was modified to a unique curve T encircling all the real solutions 
Xfc, and avoiding the complex Bethe solutions (this is possible because they are finite in 
number), as in the figure |10|. 

Clearly the T curve must be contained in the strip 

< r]+, rj^ < min{7r, np, \Im Ck\ V /c} 

Without loss of generality, assume that r]^ = r]^ = r]^ and deform T to the contour of the 
strip characterized by rj. The regions at ±oo do no contribute because of the vanishing 
of 0', then the integral can be performed on the lines /i = x ± 277, where x is real. After 



39 



X 

X 







<t 










-1 




— X — X X e 


X G e 

£=- ^ 


r 







X 



Figure 10: Contours for the integration. The crosses are roots while the circles are holes, 
algebraic manipulations involving integrations by parts and convolutions (for details see 



30|| ) one arrives at a non-linear integral equation (NLIE) to be satisfied by Zj^{€) 

Nh Ns 

cosh 9 



Z^i^) = 2N arctan — ^ + J2x{^-h)-2j2 - Vk)- 

k=l k=l 
Mc Mw 

- 52 - ^'^) - - Wk)n+ 

k=l k=l 

/oo 
dpG{d -p-ir]) log (1 + {-ife^^^^P+i^)) 
-oo 



The kernel 



1 «inh ~ 

^ ' 271 2sinh^ coshf ^ 7 

presents a singularity at the same place where 0i('(9) does. Any analytic continuation 
outside the fundamental strip < \Im'd\ < 7rmin(l,p) (I determination region) must take 
this fact into account. The source terms are given by 



X{^) = 2tt dxG{x) 
Jo 



and 



x{^) + x{^ ~ "^TTsign {Imd)) , p > 1 



is a modification of the source term due to the analytic continuation over the strip < 
|/m?9| < 7rmin(l,p), i.e. in the so called II determination region. 

Such equation is equivalent to the Bethe Ansatz, in the sense that, once solved, it gives 
the counting function from which the Bethe roots can be reconstructed. Once the Bethe 
roots are known, one can use them into eqs. ([53|) to compute the energy and momentum 
for a given state. 
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6.4 Continuum limit 



Although such NLIE is already a precious tool for the lattice model itself, its importance 
becomes essential when a continuum limit is done. 

The continuum limit must be performed carefully in the light-cone lattice setup. As 
already mentioned, one has to take N ^ oo (the so called thermodynamic limit of Sta- 
tistical mechanics) and the a — limit simultaneously, in such a way that the product 
L = Na keeps finite. This is the way to implement finite size in the continuum theory, 
removing all lattice artifacts but without destroying the cylindrical geometry. However, 
one can convince himself, by performing calculations explicitly, that if the limit is taken by 
keeping the 6 parameter fixed, the lattice NLIE blows up to infinity and looses meaning. 
This reflects the fact that the number of roots, and hence of Bethe equations goes as 
and becomes inflnite in the thermodynamic limit, thus becoming intractable. However, as 
already forecast in ref.|^, if one assumes a dependence of on A^ of the form 



4N 

e«log— . (62) 
it is possible to get a finite limit of the lattice NLIE. This limit is exactly the one that 



brings the lattice fermion introduced in |5]J into the Thirring fermion on the continuum. 
It is therefore very natural to give meaning to the continuum NLIE thus obtained as the 
remnant of Bethe Equations for the renormalized QFT. Notice that sending O ^ oo in 
this way naturally introduces a renormalized physical mass A/i. 
The continuum counting function is defined by: 

Z{^) = lim Zn{^). (63) 

— > oo 

We present here the continuum NLIE in its full generality 

/oo 
dxG{^ -x-i7]) log (1 + (_i)<5e^^(^+*'?)) (64) 
-oo 

where / = AiL. The source term g{'(}\^j) is 

Nh Ns Mc Mw 

k=l k=l k=l k=l 

The positions of {'dj} = {hj , yj , Cj , Wj} are fixed by the Bethe quantization conditions 

The parameter S can be both or 1 in principle. On the lattice it was determined by the 
total number of roots, which has become infinite here together with A^. We shall see that 
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locality implements restrictions on 5. The vacuum corresponds to a choice 6 
parameter a takes into account the twist 



0. The 



a = UJ ■ 



p+1 p — 1 



TT- 



p 



2p 



S 



+ - 



2 p+1 7r 



S UJ 



2 p+1 vr 



(65) 



In d^) the twist term uj is invariant for the shift 



UJ UJ + TX . 



The same sort of invariance is required in a for the NLIE (which is equivalent to Bethe 
equations). It is simple to verify that the expression for a ([65|) displays the following 
symmetry 

a — > a + 271 when uj —>■ uj + n . 



Note that shifting a by 27i is an invariance of the NLIE (p4\ ) , if an appropriate redefinition 
of the Bethe quantum numbers is made: Ij -^j + 1- This shift does not affect physical 
quantities, that depend only on the variables 13 j. The untwisted case corresponds to a = 0. 

In practical calculations, one starts from a Z{i)) without any roots other than the real 
ones (i.e. the counting function of the vacuum). This will be a solution of ( |6^ ) with 
g{{}\'dj) = 0. It can be got by numerical iteration starting from the driving term lsmh.{}. 
The numerical convergence is normally quite fast. It is more efficient for large /, while it 
becomes a bit slower in the region near I — > 0. Once the vacuum Z{{}) is determined, one 
can start fixing the number of roots, holes, etc... that must appear in a particular state 
he wants to study. This configuration must, of course, be compatible with the counting 
equation ([58|) . The first trial positions for G {hj , i/j , Cj , Wj} will be determined by 
the quantization conditions Z{{}j) = 27x1 j. The (integers or half-integers) numbers Ij are 
chosen in principle freely, but it turns out that some of them are constrained in terms of the 
others. In such a way a first trial of g{;d\'dj) is created and put into the NLIE. Iteration is 
redone, in order to get a better approximation for Z{'d), from which a better determination 
of the positions of can be obtained again by using the quantization rules, which in turn 
gives a better g{;d\'dj), and so on so forth, up to convergence. All this procedure can be 
done in few minutes of computing on a typical Linux/Intel platform without resorting to 
any supercomputer or other technically advanced tool. 



6.5 Energy and momentum expressions 

With a procedure analogous to the one sketched above, i.e. apply the residue trick and 
then do a continuum limit along the lines of the previous section, it is possible to give 
integral expressions for the energy and momentum. Starting from (^) one can isolate 
an extensive (i.e. proportional to A^) bulk term to be subtracted on the continuum. The 
remaining finite part of the energy can be written in a form where the positions of the 
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various kinds of roots and the function are the only inputs 

(Nh Ns Mc Mw 

cosh hj — 2 cosh yj — cosh Cj + ^^(cosh Wj)[[— 
j=i i=i j=i i=i (66) 

fir \ 
- / —2Im [sinh(x + ir/) log(l + (-l)5e'^("+'''))] 

(Nh Ns Mc Mw 

sinh hj — 2 sinh yj — sinh Cj + ^^(sinh Wj)n— 
i=i i=i i=i i=i (67) 

_ / —2Im [cosh(x + 17]) log(l + (-l)^e*^("+'''))] 
Therefore, once the Z{d) and the t?-, positions are calculated by the numerical iterative 



procedure illustrated above, it is sufficient to plug their values into the expressions ( pBD 
and (|67D to get the values of these quantities for the given state with the wanted accuracy. 
These values are numerical hut exact, in the sense that, although there is lack of a closed 
formula for them, their computation is done without any approximation other than the 
technical ones introduced by the computer truncations. 

6.6 Physical interpretation 

The limit procedure described in the previous section is mathematically consistent, but 
the question is if from the physical point of view it describes a consistent quantum theory 
and allows for a meaningful physical interpretation. 

First of all, an important remark must be made about the allowed values for the XXZ 



spin S. It is clear from (|52|) that on the lattice only integer and nonnegative values can be 
taken into account for S. But on the continuum the definition of S is no more related to 
the Bethe state (that is undefined), instead it is given implicitly by ([5^). Then, "a priori", 



there are no arguments that constrain its values to be integers. As shown in |gT|, the 
half-integer choice for S is necessary (and gives completely consistent results) to describe 
odd numbers of particles. 

At this point the following physical interpretation can be proposed. It will be refined to 
describe the correspondence with particles. Also it will be supported by many arguments 
that will be clarified in the following: 

• the physical vacuum (Hamiltonian ground state) corresponds to absence of sources 
(i.e. holes, complex roots,...); all the sources are excitations on this vacuum. Indeed, 
it can be shown that they give positive contribution to the energy, behaving as particle 
excitations on a vacuum state. 

• for = and at the various values of 5* this theory describes the sine-Gordon/massive 
Thirring model on a finite space of size L; 25* is the topological charge and can take 
nonnegative integer values. 
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• for the values 

w = — , k=l,...,q'-l (68) 

it describes the quantum restrictions of sine-Gordon model, i.e. the massive integrable 
theory obtained perturbing the minimal models Vir{r, s) by the operator $(1,3). 

• the real roots have disappeared from the counting in the continuum limit. They 
actually become infinite in number and are taken into account by the integral term, 
both in the NLIE and in the energy-momentum expression. They can be interpreted 
as a sort of Dirac sea on which holes and complex roots build particle excitations. 
Note that the presence of holes or complex roots, through the source term g^i^^l'dj) 
distorts the Dirac sea too. This can be interpreted as vacuum polarization due to 
the presence of particles. 

Observe that it has been assumed that only nonnegative values of S are required to describe 
the whole Hilbert space of the theory. Indeed the theory is assumed charge-conjugation 
invariant. Then, negative topological charge states have the same energy and momentum 
as their charge conjugate states. The assumption that all the states can be described by the 
NLIE is absolutely not trivial. A mathematical proof of this statement is not yet available, 
but a number of specific cases supports this conjecture quite strongly. 

6.7 The infrared limit of the NLIE and particle scattering 

The first task to understand the physics underlying the NLIE we have found, is to identify 
the theory as a factorized scattering, reconstructing the S-matrix. 

In the IR limit I ^ 00, the convolution term in (|6j) vanishes exponentially fast, so it 
can be dropped. Consider first a state with Nh holes only and XXZ spin S = Nh/2. 



Nh 

Z{§) = I sinh ^ + ^x{^- hj) , Z{hj) = 2TxIj (69) 

Observe that the function x can be written as 

X{^) = -i\ogSXX{d) 

where S'^^(^9) is the soliton-soliton scattering amplitude in sine-Gordon theory ||32|, if the 
parameter p = ^ — 1 is identified with the one introduced in sG theory p = gj_^2 • This 
fact is not surprising: eq.(|69D can be interpreted physically as a quantization of momenta 
in a box, for a system of particles interacting with x as logarithm of the S-matrix. The 
fact that such quantization is realized supports the interpretation of holes as solitons with 
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rapidities hj. This is further evidenced by considering the formula for the energy in such 
a case. There too the integral term drops at IR, and the formula looks like 

N„ 

E = M ^ cosh hj 
i=i 

which is the energy of Nh free particles of mass M.. The identification with the particular 
element S^^ of the S-matrix forces to give to these solitons a topological charge +1 each, 
which is consistent with the interpretation that Q = 2S. To be correct, as SZZ = S^j^ 
an analogous interpretation is possible in terms of pure antisolitons, all with charge -1, 
reflecting the charge conjugation invariance of the theory. 

Next state to consider is two holes and a complex pair. After some manipulation, the 
source terms can be arranged, thanks to some identities satisfled by x functions, in the 
form 



Z{-di) = I sinh {^i) - i log S-{^i - i^j) = 27ilj , i, j = 1,2 

where 



sinh 



= y^SXti^) 

sinh 

which is the scattering amplitude of a soliton on an antisoliton in the parity-odd channel. 
The quantum numbers /+, /_ of the two complex roots are constrained to be J± = =f| for 
consistency of the IR limit. This state has S = 0. It must have topological charge 0, again 
consistent with Q = 2S. 

It is known that there are two independent scattering amplitudes for soliton - antisoliton 
scattering (usually they are presented as transmission and reflection, here the alternative 
basis of parity even and odd channels is preferred). The even channel is realized by the 
state with two holes and a selfconjugate root 

Z{'di) = I sinh {di) - i log S+i^i - dj) = 27r/j , i , j = 1, 2 

where 



cosh 



9+iTr 



cosh 

which is the soliton-antisoliton amplitude in the parity-even channel. 

Thus, the whole soliton scattering S-matrix of sG theory has been reconstructed. Note 



that, according to the counting equation (|58D , in the repulsive regime these 3 possibilities 



exhaust the possible configurations with Nh = 2. 



45 



In the attractive regime one has also to consider the breather particles that appear as 
soliton - antisoliton bound states. It turns out that the breathers are represented by self- 
conjugate roots (1st breather) or by arrays of wide roots (higher breathers). Consistency of 
quantization at the IR limit constrains these arrays to appear only when p < 1 decreases 
under the threshold point p = 1/k for the k-th breather. Similar considerations to the 
above lead to the check that also the breather - soliton and breather - breather S-matrix 



are correctly reproduced ^ 



Nh = 1,3 states have been examined too ||4T|. In particular the one particle states 
give the correct description of a single soliton (if 5 = 1) or massive Thirring fermion (if 
6 = 0). In the former case, the particle can be at rest has quantum number Ih = 0. In the 
latter, instead, the lowest possible quantum numbers, corresponding to the lowest possible 
rapidities, are Ih = ±|. There is actually a pair of states, in accordance with the fermion 
having two components. 

6.8 UV limit and vertex operators 

This identification gets more support if we consider, for each of the states described above, 
the opposite / — > limit, where we expect to make contact with the UV limit of sG/mTh 
theory, i.e. with the c = 1 CFT described in section |3?^ . The UV calculations are usually 
more difficult to perform than the IR ones, as they involve considering the splitting of the 
NLIE at UV into two independent left and right parts (as it should be in a good CFT!), 
the so called kink equations and then expressing the energy and momentum in closed form 



in terms of these, thanks to a lemma presented in ||26|. For the details of this sort of 
manipulations the reader is invited to consult the thesis where all the calculations are 
done in detail. Here we present the main results and the physical insight they imply. 

A first important result is that the c = 1 CFT quantum number m (winding number), 
which is identified with the UV limit of the topological charge, can be related unambigu- 
ously to the XXZ spin by ±m = 2S. Of course, the ± reflects the charge conjugation 
invariance of the theory. Then, by examining the states we have already visited at IR, we 
can establish a bridge between particle states and vertex operators of c = 1 theory. 

1. The vacuum state is the one with no holes of complex roots: only the sea of real roots 
is present. There are two possible choices: 5 = or 1 (corresponding, if reinterpreted 
as a = 0,7r, to periodic and antiperiodic boundary conditions respectively). The 
result of the UV calculation gives 

for(5 = 0: A± = i.e. I 

for(5 = 1 : = i.e. V(±i/2,o) 

i.e. the physical vacuum is the one with 6 = 0. The other state belongs, with 
reference to fig. ||, to the unphysical sector IV that never contributes to make a local 
QFT at UV. 

2. The two-soliton (two-hole) state, with the minimal rapidity choice obtained with 
the lowest possible (in absolute value) quantum numbers, gives 
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(a) for (5 = and I, = -I^ = I ^ = ^ i.e. V^(o,2). 

(b) for 5 = 1 and /i = —I2 = 1 =^ a V(±i/2,2) descendent, not in UV sG spectrum, 
as it also belongs to sector IV. 

The symmetric soliton-antisoliton state (two holes and a self-conjugate root), 
always with the minimal rapidity choice 
5 = 1, Ji = -/2 = 1 and /± = ^ A± = 2^ i.e. \/(±i,o) 



The antisymmetric soliton-antisoliton state (two holes and a complex pair) 

5 = 0, Ji = /- = -I2 = -/+ = i ^ A± = 2^ i.e. V(±i,o) 
It is obvious that these last two give two linearly independent combinations of the 
operators V(±i o), one with even, the other with odd parity. 



5. The one hole state with 1 = 0,6=1 => = ^ i.e. the vertex operator V(o,i), 
belonging to sector II, thus giving more support to the 1 soliton interpretation. For 
5 = there are two minimal rapidity states with / = ±|. They are identified with 
the operators V(±i/2,±i). As these states belong to sector III, they are of fermionic 
nature and actually one identifies them with the components of the Thirring fermion. 

6. The three holes state with Ji = 1, /2 = 0, Is = — 1 =^ ^(0,3) • Higher values of 
quantum numbers give V(^n,3), n G Z 

7. The three holes and close pair state gives the first nonchiral descendent of V(o,i). 

These examples, taken all together, suggest a quantization rule (i.e a consistent choice of 6) 
to obtain states which are present in the UV spectrum of sine-Gordon or massive Thirring 
theory. It looks like 

Q + 6 + Msc e 2Z for Sine - Gordon 

S + Msc G 2Z for massive Thirring ^ ' 

where Msc is the number of selfconjugate roots. This selects sector I and II for sine-Gordon 
states, and sectors I and III for Thirring ones, as it should be from the discussion of sections 



377| and It always excludes sector IV, that does not contain local operators. All this is 
in accordance with the correct interpretation of Coleman equivalence of Sine-Gordon and 
Thirring models: even topological charge sectors are identical and the difference of the 
two models shows up only in the odd topological charge sectors, for which the content of 
Thirring must be fermionic while that of Sine-Gordon must be bosonic. 

To conclude these remarks, we briefly comment about special objects that were in- 
troduced in the classiflcation of roots but never used later. Recall their deflnition: they 
are roots or holes i/i having Z'ijji) < 0. Now, the function Z is globally monotonically 
increasing. Indeed its asymptotics for i) ±cx) are dominated by the term / sinh d which 
is obviously monotonically increasing. Also, for / large, this term dominates. Therefore at 
IR the function Z is surely monotonic and no special objects can appear. This is why we 
have not considered them in the classification of states above. However, it can happen that 
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at some critical value Icru of/, coming back from IR towards UV, the convergence of the it- 
erative procedure breaks down, thus revealing that some singularity has been encountered. 
For the scaling function to be consistently analytically continued after this singularity to 
reach the UV regime, one needs to introduce modifications to the NLIE creating exactly 
the contributions that we have called special objects. A more careful analysis reveals that 
these singularities are produced by the logarithm in the convolution term going off its 
fundamental branch. A detailed treatment of these objects can be found in [50 . 



6.9 Twisted NLIE and minimal models 



It is a well known fact ||4^ that the perturbation of the Virasoro minimal model Vir{r, s) 
by its relevant primary operator $(1,3) is integrable and is described by an RSOS restriction 
of sine-Gordon theory with 

r 

P= • (71) 

s — r 

We will use for this model the shorthand notation Vir{r,s) + $(1,3). Al. Zamolodchikov 



has put forward the idea of modifying sine-Gordon NLIE by a twist a |Q to deal with 
conformal minimal models. It is to allow this sort of approach that since the beginning we 
have considered the light-cone lattice construction with a twist u that induced a twist a 
in the NLIE. All the physics of sG/mTh model can be done with a = 0. Now we explore 
the situation where a can take well chosen nonzero values implementing the reduction and 
allowing therefore the NLIE to describe scaling functions of perturbed minimal models too. 

Looking first at the ground state, that in analogy with sine-Gordon is expected to be a 
sea of real roots, the source in NLIE is put to zero and we choose half-integer quantization 
with 6 = 0. The ultraviolet limit of the scaling functions can be computed and gives 

(72) 



p + 1 Vyr/ 

Only in the unitary models c is the Virasoro central charge. Using the rule suggested in 



(|6^) for 00 and the expression (|65D for a gives: a = vr/r. Then 



1-1 



rs 



which is exactly the effective central charge c = c — 24Amj„ of the minimal model l^ir(r, s). 
Therefore one can expect that the twisted equation describes the ground state of the model 
Vir{r,s) + $(1,3). In fact, Fioravanti et al. ||3^ calculated these scaling functions for the 
unitary case s = r + 1 and showed that they match perfectly with the TBA predictions 
already available. Moreover, choosing the following values for the twist 

a = ± — , k = l...r-l (73) 
r 
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they obtained the conformal weights of the operators $(fe,fc) , A; = 1 . . . r — 1 in the UV 
limit (the sign choice is just a matter of convention). In our notation, denotes the 

primary field with conformal weights 

The models Vir{r, s) + $(1,3) have exactly r — 1 ground states. In fact, one can see from the 
fusion rules that the matrix of the operator $(1,3) is block diagonal with exactly r — 1 blocks 
in the Hilbert space made up of states with the same left and right primary weights. In 
each of these blocks, there is exactly one ground state and for the unitary series s = r + 1, 



it was conjectured in |^ that their UV limits are the states corresponding to One 
can check that in the general non-unitary case the twists ([73D correspond in the UV limit 
to the lowest dimension operators among each of the r — 1 different blocks of primaries 
(see explicit examples later). 

These ground states are degenerate in infinite volume, but for finite / they split; their 
gaps decay exponentially as / — *• 00. In the unitary case, they were first analyzed in the 
context of the NLIE in where it was shown that the NLIE predictions perfectly match 
with the TBA results already available for the unitary series. 



Ground states for non-unitary models have been treated in [0. The simplest example 
is the scaling Lee- Yang model Vir{2, 5) + $(1,3), for which 

Mb = 2M sin ^ = VSM 

is the mass of the fundamental particle of the Lee- Yang model (this is more natural here 
than using the mass M. of the soliton of the unrestricted sine-Gordon model as a scale, 
since the soliton disappears entirely from the spectrum after RSOS restriction). Call 

Ib = MbL = V3l 

where / is the variable appearing in NLIE. There is only one independent value of the twist 

TT 

There is only one ground state in this model, which corresponds to the primary field 
with conformal weights 

A = A = -1 , 
5 

All the models of the class Vir{2, 2n + 1) + $(1,3) have only one ground state. For models 
with two ground states, we can take a look at Vir{3,5) {Vir{3,7) was also taken into 
account in ||42|). For Vir{3,5) the ultraviolet spectrum is defined by the following Kac 
table, where the weight of the field ^(k,i) is found in the k-th row and l-th column. 
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-1/20 


1/5 


3/4 


3/4 


1/5 


-1/20 






The two blocks of the perturbing operator $(1,3) are defined by the fields {$(1,2) , "^(1,4)} 
and {$(1,1) , $(1,3)}, respectively. The ground states correspond in the UV to the operators 
$(1,2) and $(1,1), as can be checked directly using formulae (|72D , (|7lD and ([73|). 

Now we turn our attention to excited states over these minimal model vacua. We 
restrict ourselves to the case of neutral (i.e. S = 0) states. Even for states with a zero 
charge the relation between a and uj is highly nontrivial. 

It is quite easy to show that choosing the value of u as 



p+l 



where k is integer, we can reproduce all the required values of a listed in equation ([73|). 



The twisted lattice Bethe Ansatz was analyzed by de Vega and Giacomini in On the 
lattice, passing from the sine-Gordon model to the perturbed Virasoro model amounts to 
going from the six-vertex model to a lattice RSOS model. In it was shown that to 



obtain all the states of the RSOS model it is necessary to take all the twists 

kir 



LU 



P+l 



mod TT 



into account. The fact that not all these twists correspond to inequivalent values of a and 
so to different physical states is a consequence of the RSOS truncation. 

We remark that the parameter a drops out of the second determination of Z in the 
attractive regime. This is important because as a consequence the IR asymptotics of the 
breather states do not depend on a and so the S-matrices involving breathers are un- 
changed. In fact, scattering amplitudes between solitons and breathers remain unchanged 
too. This matches with the fact that the RSOS restriction from sine-Gordon theory to per- 
turbed minimal models does not modify scattering amplitudes that involve two breathers 
or a breather and soliton 
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Let us start by examining the scaling Lee- Yang model Vir{2, 5) + $(1,3). Since we have 
a single ground state, there are no kinks in the spectrum. We fix the value of a as above, so 
we still have a freedom of choosing 6. This can be done by matching to the UV dimensions: 
if for a certain state we choose the wrong value of 6, we find a conformal dimension that 
is not present in the Kac table of the model. 

The excited states are multi-particle states of the first breather of the corresponding 
unrestricted sine-Gordon model, which has P = §• One can calculate the state containing 
one particle at rest. It turns out that as we decrease /, the self-conjugate root starts 
moving to the right. It does not remain in the middle like in the a = case, which is 
to be expected since for nonzero a we have no left/right symmetry. However, the total 
momentum of the state still remains zero due to a contribution from the integral term in 
momentum equation ( |6^ ) . We have here an example of the phenomenon of the appearance 
of the special root and its two accompanying holes: the numerical iteration breaks down 
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at around Icrit = 2.5. For / < Icru the NLIE must be corrected by the introduction of this 
special root and the two accompanying holes (which are not special). They must be taken 
into account to compute the UV limit correctly. Notice that the two holes thus introduced 
do not have any counterpart at IR and it is not possible to give them any interpretation as 
particles. Their quantum numbers cannot be freely chosen. They are constrained to take 
the same value of the quantum number of the special root, which in turn is fixed by the 
constraints of the problem. 

We know that the Lee- Yang model contains only two primary fields, the identity I and 
the field ip with left /right conformal weights -1/5. In fact, the ground state of the massive 
model corresponds to if in the UV limit. For the one particle state, it turns out that the 
special root and one of the holes move to the left together with the self-conjugate root, 
while the other hole moves to the right. This gives, after the "kink" NLIE calculations are 
done, A = A = 0, i.e. the identity operator I, as expected from TCSA in |3^ . 

Let us look now at moving breathers. If the self-conjugate root has Bethe quantum 
number / = 1, the corresponding state will have momentum quantum number 1, i.e. 

p = H^ 

R • 

and in the UV A — A = 1. No special root appears here. The reason is that the self- 
conjugate root moves to the left and the real part of its position ^ is given to leading order 
by 

27r/ 



sinh(3f?e^?) 



Ib 

As a result, the contribution to the derivative of Z from the Zsinht? term remains finite 
when / ^ 0. In the previous example of the particle at rest the left-moving nature of the 
self-conjugate root when 1 = does not prevent the occurrence of the breakdown in the 
iteration scheme: since its Bethe quantum number is zero, it does not move fast enough 
to the left in order to balance the negative contribution to derivative of Z coming from 
the self-conjugate root source. At the moment we have no way of predicting analytically 
whether or not there will be specials in the UV limit: we just use the numerical results to 
establish the configuration for the evaluation of UV weights, supplemented with a study 
of the self-consistency of the solution of the kink equation. The UV dimensions for the 
moving breather turn out to correspond to the state L^iip. 

One can similarly compute the UV dimensions for some other excited states. For 
example, the two-particle states with half-integer Bethe quantum numbers Ji > 0, /2 < 
for the two self-conjugate roots are found to have 

A = h/i + -, A = /2 + -, 

5 2 ' 5 2 ' 

which shows that they correspond in the UV to descendent states of ip. The first such state 
with quantum numbers 

2 ' 2 
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corresponds to L^iL^iip. 

The lowest lying three-particle state of zero momentum, with Bethe quantum numbers 
(—1,0, 1) corresponds to the left/right symmetric second descendent of the identity field, 
i.e. to the field TT, where T denotes the energy-momentum tensor. This is very interesting, 
since from experience with NLIE UV calculations one would naively expect this to be a 
first descendent (descendent numbers are usually linked to the sum of Bethe quantum 
numbers of left/right moving particles and this state is the lowest possible descendent of 
the identity I). However, the field L_iL_iI is well-known to be a null field in any conformal 
field theory. It seems from this example that NLIE is clever enough to avoid null vectors 
in minimal models, but of course more investigation is needed in this respect. 

The above correspondences are confirmed by comparing to TCSA data (see the won- 
derful figures in In general, one can establish the rule that states with odd number 
of particles must be quantized by integers {6 = 1), while those containing even number of 
particles must be quantized by half-integers (5 = 0) in order to reproduce correctly the 
spectrum of the scaling Lee- Yang model. 

We conducted similar studies for the models Vir{2, 7) + $(1,3) and Vir{2, 9) + $(1,3) and 
found similarly good agreement with TCSA data. For the first one-particle state of the 
model Vir(2, 7) + $(1,3) we also checked our results against the TBA data in the numerical 
tables of [3^ and found agreement with the TBA results. 

Given the choice of a above, the correct rule of quantization in all of the models 
rir(2,2n + l) + <l>(i,3) is 

S + Msc^ 2Z 



where Mgc is the number of self-conjugate roots in the source corresponding to the state. 



This is exactly the same rule as the one established for pure sine-Gordon theory in ||28[. In 
the presence of the twist, such a rule of course has meaning only together with a definite 
convention for the choice of a. 

It is interesting to note that in the case of Vir{3, n) models, all the neutral states must 
come in two copies, since they can be built on top of either of the two ground states. We 

take the example of the Vir{?), 7) model and the states corresponding to a breather at rest. 

3 

4' 



We have p = |, but now there are two inequivalent values for the twist 



TT 2tI 

" " 3 ' T' 

When a = 0, we can estimate the critical value of / where specials appear to be lent = 5.23. 
In this case, the twist helps a bit, because it makes the self-conjugate root a left mover; it 
is intuitively clear that the bigger the twist, the more it lowers the eventual value of Icrit- 
From TCSA data one can identify that breather #1 is really the one-particle state in the 
sector of ground state #1 (a = 7r/3), while breather #2 is in the sector over ground state 
#2 (a = 27r/3). 

A direct calculation of the conformal weights gives ^ for breather #1 and for breather 
#2, which are in complete agreement with TCSA data. We also checked the two different 
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states containing two breathers with Bethe quantum numbers Ji = ^ , J2 = — ^ , and 
found an equally excellent numerical agreement with TCSA. Just like in the case of sine- 
Gordon and scaling Lee- Yang model, for these states one can continue the iteration of 
the NLIE down to any small value of /, although at the expense of a growing number of 
necessary iterations to achieve the prescribed precision. 



7 Conclusions 

We have shown how the nonlinear integral equation deduced from the light cone lattice 



model of ||22|| describes exactly the Finite Size Effects of the sine-Gordon/massive Thirring 
theory, taking into account excited states too. The most important results are summarized 
as follows: 

1. By examining the infrared limits of the equation it has been shown that it leads to 
the correct two-particle S-matrices for both scattering states and bound states; 

2. By computing the UV conformal weights from the NLIE we have shown that the 
energy-momentum spectrum is consistent with the UV one of sG/mTh theory only 
if we choose the parameter 6 (i.e. the quantization rule) as indicated in ([70|). 

3. The predictions of the NLIE have been verified by comparing them to results coming 
from the TCSA approach (for sG/mTh). 

4. The framework required to deal with minimal models perturbed by $(1,3) is built up. 
Many examples and numerical/analytical checks are given for the ground state in the 
unitary p8[ and non unitary cases. 



All the conformal dimensions can be reproduced (Kac table) with a convenient choice 



of the twist (|73|) . The particular relation suggested in is not enough to describe 



the whole spectrum. 

6. The IR computations reproduce correctly the S-matrix of perturbed minimal models 
(at list in the cases not involving kinks), because the p < 1 second determination 
drops the twist. 

7. Numerical calculations (TCSA) of concrete examples and comparison with TBA 
where possible, give a strong evidence for the correctness of the energy levels de- 
rived from the twisted NLIE for excited states, both for the full sG/mTh model and 
for the Restricted models. 

The understanding of "(twisted) sine-Gordon NLIE" and of the finite size behaviour of 
the continuum theory defined from the NLIE (|64|) is not quite complete. Indeed some 
open questions have not an answer yet and further interesting developments may be fore- 
cast. We close this review by pointing to possible interesting questions that need further 
investigation. 
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1. Is the set of scaling functions provided by the NLIE complete i.e. can we find to every 
sG/mTh or minimal model state a solution of the NLIE describing its finite volume 
behaviour? The main difficulty in this "counting problem" is that the structure of 
the solutions is highly dependent on the value of the coupling constant - consider 
e.g. the appearance of special sources. 

2. The multi-kink states characteristic of the perturbed minimal models (except for the 
series Vir{2, 2n + 1)) have been omitted in the analysis performed in the last section. 
Although the general treatment is valid for such states too, a detailed description is 
far more complicated than for states which contain only breathers and is left open 
to further studies. 

3. There is also an unresolved technical difficulty, namely that the source configuration 
of the NLIE may change as we vary the volume parameter /. Typically what hap- 
pens is that while the counting function Z is monotonic on the real axis for large 
volume, this may change as we lower the value of / and so-called special sources 
(and accompanying holes) may appear. We do not as yet have any consistent and 
tractable numerical iteration scheme to handle this situation, although the analytic 
UV calculations and intuitive arguments show that the appearance of these terms in 
the NLIE is consistent with all expectations coming from the known properties of 
perturbed CFT. In addition, in the range of / where we can iterate the NLIE with- 
out difficulty, our numerical results show perfect agreement with TCSA. We want 
to emphasize that these transitions are not physical: the counting function Z and 
the energy of the state are expected to vary analytically with the volume. It is just 
their description by the NLIE with requires a modification of source terms. As it was 



pointed out also in ||23, |2^, the whole issue is related to the choice of the branch of 



the logarithmic term in the NLIE (|64|). 



4. The problem pointed out at the previous point is very similar to the behaviour of sin- 
gularities encountered in the study of the analytic continuation of the TBA equation 
[p^ and we can hope that establishing a closer link between the two approaches can 
help to clarify the situation. From the form of the source terms in the NLIE it seems 
likely that the excited state equations can be obtained by an analytic continuation 



procedure analogous to the one used in TBA [3§| to obtain the excited state TBA 
equations. Certain features of the arrangement of the complex roots in the attractive 
regime and their behaviour at breather thresholds also point into this direction. This 
is an interesting question to investigate, because it can shed light on the organization 
of the space of states and can lead closer to solving the counting problem described 
above. 

5. Even if on the lattice the Bethe vectors given in (|^) completely describe the Hilbert 
space of the XXZ chain, the form of eigenvectors for continuum energy and momen- 
tum is completely unknown. This fact refiects the unresolved question of the form 
of energy and momentum eigenvectors for the minkowskian sine-Gordon theory (the 
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Faddeev-Zamolodchikov algebra is only a phenomenological picture). Also the deter- 
mination of correlation functions is in an early stage, although in recent publications 



IH| explicit expressions for the field (e*""^) and its descendent are given. 



6. The extension of light-cone approach and NLIE to other QFT models is in progress. 
An interesting extension is the description of finite volume spectrum of Zhiber- 
Mikhailov-Shabat model (imaginary Bullough-Dodd). A first suggestion in this di- 
rection, even if is obtained in a completely different approach, recently appeared in 



the [^, for the ground state. Early attempts to get the excited states can be found 
in Extensions to vertex models based on higher rank simply-laced Lie Algebras 



have also been investigated |51, should correspond - a detailed analysis has 



not been done yet - to the solitonic (imaginary coupling) affine Toda Field Theories. 

The second and fourth points are important because their investigation may lead closer to 
understanding the relation between the TBA and the NLIE approaches. It is quite likely 
that establishing a connection between the two methods would facilitate the development 
of both and may point to some common underlying structure. 

As a final remark, let us mention an application of the NLIE to a chemical compound, 
i.e. "copper Benzoate" (Cu(C6D5COO)2-3D20). Its specific heat has been computed in ||19| 
using the vacuum untwisted (a = 0) NLIE, where the finite size L is the inverse temperature 
L = {kBT)~^ . This corresponds to do equilibrium thermodynamics of sine-Gordon model. 
It is curious to notice that in that context the NLIE is called Thermal Bethe Ansatz, 
emphasizing the well known relation between statistical mechanics and QFT. Exactly in a 
similar contest (Heisenberg ferromagnet model) the original Bethe Ansatz |^ was born. 
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